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Geometric techniques have played an important role in the seventies, for the study of the spectrum 
of many-body Schrodinger operators. In this paper we provide a formalism which also allows to 
study nonlinear systems. 

We start by denning a weak topology on many-body states, which appropriately describes the 
physical behavior of the system in the case of lack of compactness, that is when some particles are 
lost at infinity. We provide several important properties of this topology and use them to write 
a simple proof of the famous HVZ theorem in the repulsive case. In a second step we recall the 
CO ' method of geometric localization in Fock space as proposed by Derezihski and Gerard, and we 

^ ' relate this tool to our weak topology. 

, We then provide several applications. We start by studying the so-called finite-rank approxima- 

CO i tion which consists in imposing that the many-body wavefunction can be expanded using finitely 

many one-body functions. We thereby emphasize geometric properties of Hartree-Fock states and 
C^l ' prove nonlinear versions of the HVZ theorem, in the spirit of works of Friesecke. 

In the last section we study translation-invariant many-body systems comprising a nonlinear 
term, which effectively describes the interactions with a second system. As an example, we prove 
the existence of the multi-polaron in the Pekar-Tomasevich approximation, for certain values of 
the coupling constant. 
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Introduction 

A system of N (spinless) quantum particles is usually described by an energy functional 
^ H> £(^) G K where ^ is a normalized function of the A-body space 

N 

?) N ■■= L 2 {R d ) ~ L 2 ((R d ) N ). (0.1) 

Here d is the dimension of the space in which the N particles evolve, that is d — 3 in the 
physical case. If the particles are indistinguishable bosons (resp. fermions) it is additionally 
assumed that ^ is symmetric (resp. antisymmetric) with respect to exchanges of variables 

(x u ...,x N ) e (R d ) N . 

In the simplest case the energy £ is the quadratic form associated with a self-adjoint 
operator on Sj N . For nonrelativistic particles interacting with a two-body potential W and 
submitted to an external potential V, the corresponding N-body Hamiltonian reads 

N ( A \ 
H V (N)=J2 \-^ L + V (xj) + E W(x k - Xi ). (0.2) 

j=l ^ ' l^k<i^N 

The study of the properties of self-adjoint operators of this form has a long history [26] 
and it is certainly one of the most significant successes of mathematical physics in the past 
decades. Of particular interest is the spectrum of H V (N). 

The advent of geometric methods in the late seventies has been particularly important. 
By 'geometric' it is usually meant the use of clever partitions of unity in configuration space 
in order to relate local properties of H v (N) (seen as a partial differential operator) and 
spectral properties. Initiated in the sixties by Zhislin [66] and Jorgens and Weidmann [27] . 
the systematic use of geometric ideas in Schrddinger operators theory really started in 1977 
with the works of Enss [16] ■ Deift and Simon [11], and Simon [57]. It was then further 
developed by Morgan [45] , Morgan and Simon [46] , and Sigal [54] [55] [56] . For a review of 
these techniques we refer for instance to [52, 19] [26]. 
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A famous example of the use of geometric methods is the so-called HVZ Theorem of 
Zhislin |66j , Van Winter [62] and Hunziker |25) . Under suitable decay assumptions on V and 
W, it relates the bottom of the essential spectrum of H V (N) to the ground state energy of 
systems with less particles: 

inf a css (H v (N)) = inf {E v (N - k) + E°(k), jfe = 1, JV} , (0.3) 

where E v (N) := inf a(H v (N)) is the ground state energy for N particles. Physically this 
result says that in order to reach the bottom of the essential spectrum one has to remove k 
particles from the system and place them at infinity. The total energy is then the sum of the 
ground state energy E (N — k) of the N — k remaining particles plus the energy E°(k) of 
the k particles at infinity. The number k of particles to extract is chosen such as to minimize 
the total energy obtained by this procedure. A consequence of (|0.3|) is that E v (N) is an 
isolated eigenvalue if and only if 

E V (N) < E V (N -k) + E°(k), Vfc = 1, ...,iV. (0.4) 

Although physically quite natural, the HVZ formula (|0 . 3[) is mathematically not obvious, 
in particular because the three problems corresponding to having N, k and N — k particles 
are posed on the different Hilbert spaces S) N , S) k and ¥) N ~ k . When proving (|0.3p . geometric 
methods indeed make a crucial use of the fact that the many-body space has the structure 
of a tensor product, that is Sj N ~ $) N ~ k ® S) k . 

Linear problems are not the only possible ones occurring in the study of many-body 
quantum systems. Indeed, most numerical methods used by physicists and chemists resort 
to nonlinear models. Sometimes the energy is kept linear but the set of states is reduced by 
assuming that the wavefunctions belong to a well-chosen manifold. In some other cases 
it is convenient to modify the many-body energy £ by adding nonlinear empirical terms in 
order to account for involved physical effects which are too complicated to describe in a 
precise manner. 

Nonlinear methods also have a long history, in particular within the field of partial 
differential equations. Loosely speaking, a typical question is to understand the behavior of 
sequences of functions {ip n } (say in L 2 (R d )), in particular in the case of lack of compactness, 
that is when ip n — 1 tp weakly in L 2 but (p strongly. The sequence {ip n } can be a 

minimizing sequence of some variational problem or a Palais-Smale sequence |61j (in these 
cases the goal is often to prove by contradiction that it must converge strongly). Or it can 
be the solution of a time-dependent equation, which experiments a dispersive or a blow-up 
behavior in finite or infinite time (in this case lack of compactness has some physical reality) . 

The first to tackle such issues on a specific example were Sacks and Uhlenbeck [53] in 
1981 who dealt with a concentration phenomenon for harmonic maps. Brezis and Niren- 
berg [8] then faced similar difficulties for some elliptic partial differential equations with 
a critical Sobolev exponent. In 1983, Lieb proved in [3 3) a useful lemma dealing with 
lack of compactness due to translations in the locally compact case. A general method 
for dealing with locally compact problems was published by Lions (38j [39] in 1984 under the 
name "concentration-compactness". Later in 1984-85, Struwe [60] and, independently, Brezis 
and Coron [7] have provided the first "bubble decompositions', whereas Lions adapted his 
concentration-compactness method to the nonlocal case [40] [41] . For a review of all these 
techniques, we refer for instance to [61] . 
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When studying the compactness of minimizing sequences for a variational problem of 
the general form 

I(N) = inf £(<p), 

a useful argument is to rely on so-called binding inequalities 

I(N)<I(N-X)+I°(X), V0<\^N, (0.5) 

where I°(N) is the ground state energy when the system is sent to infinity (that is when 
all the local terms have been dropped in the energy £). Imagine that one can prove that 
a non-compact minimizing sequence {<p n } would necessarily split into pieces in such a way 
that the total energy becomes the sum of the energies of these pieces. Then an energetic 
inequality like (|0.5p yields a contradiction and implies that all minimizing sequences must 
be compact. Arguments of this type are ubiquitous in studies of nonlinear minimization 
problems. 

The formal link between the HVZ formula (|0 . 3[) and binding inequalities of the form 
of (|0.5[) has been known for a long time. There are important differences, however. In the 
HVZ case one has a quantized inequality (|0.4[) in which only an integer number of particles 
can escape to infinity. On the contrary the binding inequality (|0.5p is not quantized since in 
L 2 (R d ) the sequence {fn} can split in pieces having an arbitrary mass. Vaguely speaking, 
this comes from the fact that in the case of lack of compactness, ip n usually behaves as a 
sum of functions whereas an TV-body wavefunction is rather a tensor product. 

The goal of this paper is to present a theory which combines nonlinear and geometric 
techniques, with the purpose to study some many-body systems involving nonlinear effects. 
A first attempt in this direction was already made by Friesecke in his paper 19 on multi- 
configuration methods, a work which partly inspired the present paper. However, instead of 
concentrating only on some specific examples, a large part of this article (Sections [2] and is 
devoted to the presentation of a simple but general theory which, we hope, will be reusable 
in many other situations. We apply it to some nonlinear models in Sections 2] and [5] 

In this work, we are particularly interested in finding an appropriate description of the 
possible lack of compactness of many-body wave functions. As we now explain, usual methods 
of nonlinear analysis are rather inefficient in this respect. Consider for instance a sequence 
of two-body wavefunctions of the form: 

y n = l p®<p n , (0.6) 

that is ty n (xi,X2) — (p(xi)ip n (x2), with <p,<p n £ L 2 (R d ). We assume that ip n — * weakly in 
L 2 (Mr), hence we may think of as describing a system of two particles, one in the fixed 
state ip and the other one 'escaping to infinity'. It is then easily verified that 

* n weakly in L 2 (R d ) ® L 2 (R d ) ~ L 2 ((M d ) 2 ), 

which suggests that looking at weak limits of two-body wavefunctions does not say much on 
the real behavior of the system. We would rather like to have, for obvious physical reasons, 
that 

"*„ ^ (0.7) 

since one particle is lost and the other one stays in the one-particle state ip. However this 
does not make much sense as such, since \I/ n £ L 2 (R d ) <£> L 2 (R d ) and <p £ L 2 (R d ) live in 
different Hilbert spaces. 
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In Section [2] we introduce a very natural topology on many-body states, which we call 
geometric topology, and for which (|0 . 7[) is actually correct. The geometric topology is very 
different from the usual weak topology (as can already be seen from the fact that $ n — 
weakly). It is however the one which is physically relevant for many-body systems. 

Let us vaguely explain how the geometric topology is defined. As is suggested by (10. 7|) . 
even if we start with a sequence of states containing A particles (in the A-body space Sj N ), 
we have to allow limits in spaces with less particles. All the particles could even be lost in 
the studied process, in which case we would end up with the vacuum. For this reason, the 
behavior of A-body states must be studied in the so-called truncated Fock space 

F< N -^CeSj 1 ®---®^ (0.8) 

which gathers all the spaces of k particles, with ^ k ^ A. As we shall see on specific 
examples, it is also natural to allow a geometric limit which is a mixed state, even when 
the sequence is only made of pure states. Let us recall that a mixed state V on F^ N is 
a trace-class self-adjoint operator such that T ^ and Tr^^F) = 1. A pure state is a 
rank-one projector, r = |<&)(<&| with $ € F^ N (for instance $ = 0©---©0ffl , I'in the case 
of a pure A-body state 'F £ Sj ). 

The geometric topology on mixed states on F^ N is defined by means of the weak topolo- 
gies of all the corresponding density matrices, which are specific marginals (partial traces) 
reflecting the tensor product structure of the ambient Hilbert space (hence the name 'geo- 
metric'). The definition of the density matrices is recalled in Section [T] below. In particular 
we say that T n —^ g T geometrically when all the density matrices of T„ converge to that 
of T, weakly-* in the trace class. Let us emphasize that the geometric limit F is always a 
state, that is it satisfies Tr(L) = 1. There is never any loss in the trace norm when passing 
to geometric limits. 

For instance the one-body density matrix of our two-body sequence {'Fn} m f|0 - 6[) is the 
operator acting on L 2 (M. d ) 

(we assume for simplicity that (p _L (p n for all n). By the weak convergence of ip n 0, it 
holds 

*n ~~ IvXvl weakly-*. 

The operator | I is precisely the one-body density matrix of the one-body state tp g $y. 
We indeed have that 

© © \<p<%> <p n ){tp® <f n \ - s ' © \ip)((fi\ © geometrically in F^ 2 , 

which is the precise mathematical meaning that we can give to (|0.7p . 

Our weak topology is the restriction to states on F^ N of a well-known weak-* topology 
associated with the CAR/CCR algebra (Remark 12 .lj) . But, to our knowledge, the usefulness 
of this notion of convergence for many-body problems has never been pointed out in the 
literature. As will be seen on several examples in this work, it is however the most natural 
weak topology for many-body states. It is a crucial notion when strong convergence does 
not hold a priori, that is in the case of possible lack of compactness. 

In Section 12.1.21 we proceed to give important properties of geometric convergence. We 
start by showing that the set of states is compact for the geometric topology in Lemma l2~2l 
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This means that any sequence of states {r„} on the truncated Fock space has a 

subsequence such that T Uk ^ g T geometrically. This result is very important in applications. 
We then show in Lemma 12.31 that strong convergence is equivalent to the conservation of 
the total average particle number. 

We illustrate the use of our theory in Section 12.21 We consider an TV-body Hamiltonian 
of the form of (|0.2j) with W ^ and we show that, in contrast with the usual weak topology 
of Sj N , the associated quantum energy is lower semi-continuous for the geometric topology. 
This enables us to provide a very simple proof of the HVZ Theorem, in this particular case. 

Equipped with a new weak topology, we then need a second important notion: geometric 
localization. As we have already mentionned, localization has always played an important 
role in the study of Schrodinger operators. As we want to find out where are the particles 
which stay and where do go those which escape to infinity, we need to be able to describe 
the state of our system in a given domain D C K d . 

If we think of a one-body state ip £ L 2 (R d ), then the corresponding localized state in 
a domain D should clearly be described by the function \d<P- However, \n<p is in general 
not a state since J D \(p\ 2 < 1, except when ip has its support in D. Having removed what 
is outside D corresponds, in our language, to the vacuum state. Thus the localized state 
should rather be 



in the truncated Fock space J 7 ^ 1 . 

The correct notion of localization of any mixed state of J 7 '^ N which generalizes (|0 .9[) was 
introduced by Derezinski and Gerard [13] in the context of Quantum Field Theory. It is 
even possible to define a localization with respect to any operator B on i 2 (R d ) such that 
BB* $C 1, not only for B the multiplication operator by the characteristic function trj (this 
is in particular useful when dealing with smooth cut-off functions). In Section|3]we recall the 
definition of geometric localization in our context and we provide several of its properties. Of 
particular interest is the fact that if T n V geometrically and {r„} has a bounded kinetic 
energy, then one gets a strong convergence of the localized states in any bounded domain D. 
This generalizes the well-known Rellich compactness embedding theorem in Sobolev spaces, 
to the setting of many-body states and geometric topology. 

In Section [321 using both geometric convergence and localization we are able to provide 
a simple proof of the HVZ theorem in the general setting (when W has no particular sign) , 
which particularly enlightens a crucial but simple geometric property of V-body functions, 
see Equation (|3.6|) below. It is not our intention to pretend that our proof of the HVZ 
theorem is better than any of the other existing proofs. We rather aim at accustoming the 
reader to the techniques that we will use for nonlinear models in Sections [4] and and for 
which usual linear methods are inappropriate. 

We turn to the study of nonlinear models in Sections 0] and [5] 

In Section [4] we study the so-called finite rank approximation in which one restricts to 
V-body states which can be expanded using a finite number of (unknown) one-body orbitals 
ipx, (p r . In the bosonic case we obtain the Hartree model for r — 1. In the fermionic case 
the Hartree-Fock theory [35] is obtained when r = N (the number of particles) whereas 
r > N leads to multiconfiguration methods 19, 30 . Despite the fact that these methods are 
essentials tools of quantum physics and chemistry, their geometric properties have deserved 
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little interest in the literature so far. In [19], Friesecke was, to our knowledge, the first to 
consider both the Hartree-Fock and the multiconfiguration theories as real iV-body models 
and to use geometric techniques in order to derive nonlinear HVZ-type results. 

Our goal is to emphasize geometric properties of finite-rank states, that is to find what 
can be said on geometric limit points of sequences or on geometric localization of such 
special states. For instance we show in Section [4.1.21 that the geometric limit of a sequence 
of pure Hartree-Fock states is always a convex combination of pure Hartree-Fock states, see 
Example l4.3l below. Using such properties, we are able to provide a simple proof of Friesecke's 
results, as well as to derive other theorems. For instance in Theorem 14.21 below, we prove a 
nonlinear HVZ-type result for a translation-invariant Hartree-Fock theory, combining ideas 
of Lions {38} [39! an< ^ geometric techniques. This result is in the same spirit as what was done 
for neutron stars in a recent collaboration with Lenzmann [29 . 

In Section [5] we study another kind of nonlinear models where all possible many-body 
states are considered but nonlinear effective terms are added to the quantum energy £ 
in order to describe some specific physical effects. To be more precise we concentrate on 
translation-invariant models of the form 



where F is a concave nonlinear function of the charge density py, and H (N) is the TV-body 
Hamiltonian (|0.2[) with V = 0. In practice the purpose of the nonlinear term F(p^) is to 
model the interaction of the N particles with a second complicated system. For instance we 
consider in Section [5.2l the multi-polaron in the Pekar-Tomasevich approximation. This is a 
system of N nonrelativistic electrons with an effective nonlinear term 



modeling interactions with the phonons of a polar crystal in the regime of strong coupling. 
We show the existence of bound states for all a > t c (N) where t c (N) < 1, which covers the 
physical case. This complements recent results of jTS] [T7] . 

The paper is organized as follows. In Section [T] we provide necessary notations and 
some preliminary results. The reader at ease with the concepts of Fock space, creation and 
annihilation operators and density matrices may want to skip most of the material of Section 
[TJ Of importance is Lemma fl . 1 1 which provides crucial properties of density matrices for states 
on a truncated Fock space J 7 ^^. Section [2] is devoted to the definition and the derivation of 
important properties of the geometric topology and convergence. This is followed by a proof 
of the HVZ theorem in the repulsive case. In Section [3] geometric localization is defined and 
the general HVZ theorem is proved. Sections |4] and [5] are respectively devoted to the study 
of the finite-rank approximation, and of nonlinear systems of the form (|0.10|) . 

For the sake of clarity we usually do not state the most general results and rather favor 
some chosen applications. Many of our theorems can be generalized in several directions. 

Acknowledgment. I started this work after several interesting discussions with Enno Lenzmann. 
He was the first to draw my attention to the multi-polaron model, which was the starting point of 
this article. I am also indebted to Vladimir Georgescu for many stimulating discussions. 

The research leading to these results has received funding from the European Research Council 
under the European Community's Seventh Framework Programme (FP7/2007-2013 Grant Agree- 
ment MNIQS no. 258023). 
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1. Notation and preliminaries 

We start by fixing some important notation and vocabulary, as well as by providing some 
preliminary results that will be useful throughout the paper. The reader acquainted with 
Fock spaces can jump to Section 11.41 where density matrices arc dchncd and some of their 
important properties derived. 

1.1. Spaces and algebras 

For a (separable) Hilbcrt space F), we denote by B(S)) and JC(S)) the algebras of, respectively, 
bounded and compact operators on Sj. The Schatten space & p (f)) C K,(Sj) is defined [58] by 
requiring that 

U\\e P m :=Tr(|AH 1/p < oo, 

with \A\ = \[A*A. Operators in have a well-defined trace Tr(A) = ^\ (fi,Afi) (for 

any orthonormal basis {/,} of Sj). Operators in &2($)) are called Hilbert-Schmidt. We recall 
that [58] 

= Si(i5) and (6i(j5))'=B(i5). (1.1) 

Since K,(Sy) is separable when f) is separable, (|1.1[) means that any bounded sequence {T ra } in 
©i(fj) has a subsequence which converges weakly-* in the sense that lim.n_j.oo TT^(T n K) = 
Trfi(TK) for all K £ !€(£)). The same holds for bounded sequences in B($j), with IC(S)) 
replaced by 

In the whole paper we fix as space for one quantum particle io = L 2 (R d ). We could as 
well work in a domain Q with appropriate boundary conditions, use a discrete model, or 
even, for most of our results, take an abstract Hilbert space. These obvious generalizations 
are left to the reader for shortness. 

Similarly, for simplicity we almost always restrict ourselves to the case of quantum sys- 
tems made of one kind of indistinguishable particles (fermions or bosons) without spin. Most 
results can be easily generalized to the case of several kinds of particles having internal de- 
grees of freedom. The space for N indistinguishable fermions is the antisymmetric tensor 
product 

Si? :=f\Sj = Ll((R d ) N ) 
i 

consisting of wavefunctions ^ which are antisymmetric with respect to exchanges of vari- 
ables: V_(_i, Xi, Xj, xn) — — \-(xi, Xj, Xi, xjsr), with Xk £ K d for k = 1, N. 
The space for TV indistinguishable bosons is the symmetric tensor product 

N 

:=\/^ = ^((KT) 
1 

consisting of wavefunctions \_ which are symmetric with respect to exchanges of variables: 
...,Xi,...,Xj,...,XN) = ^(xi, ...,Xj, ...,Xi, ...,x N ), with x k £ K d for k — 1, N. 
The corresponding fermionic or bosonic Fock space is denoted as 

N>1 
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Saying differently, it is the space composed of sequences of the form ^ = (tfj° , ip 1 , ip 2 , ...) £ 
C x f) x Sy^ . x • • • satisfying the constraint that 



It is a Hilbert space when endowed with the scalar product 



The vacuum state is by convention defined as 17 := (1,0, 0, ...) G ^"a/s- 

As we consider iV-body systems, we most always work in the 'truncated' Fock space 



N 



/a 

n=X 



(1.2) 



which we identify to a closed subspace of J~ a / S . Similarly, any iV-body vector of fy^/s can be 
viewed as a vector of J-^/s or °^ 3~ a / s . As we explain later in Section [2l the 'geometric' limit 
of a sequence (ip n ) C wn ^ always live in the truncated Fock space 

For ipi € i}^ 1 and ^2 € -Q^ 2 , we define the antisymmetric tensor product tpi A V^2 € 
Sj^ +N2 as follows: 

■0! A 1p 2 {%l, -,XNi+N 2 ) 



y/mi n 2 \{N! + N 2 y. 

x ^2 s g n ( cr )V'i(^(i),---,a; C 7(Af 1 )) '02(a;<T(w 1 +i), •■•,a; C 7(jv 1 +JV 2 ))- (1.3) 

Here <S/v is the group of permutations of {1, N}. When {/,} is an orthonormal basis of 
ft, then {f h A ••• A j^}^ 

< - <ijv forms an orthonormal basis of $j a . 
For bosons, we define similarly, for ipi G -fi^ 1 and V>2 € fj^ 2 , 

?/>i V^ 2 (a:i,..-,a;jVi+JVa) := 



x /Afi!A^ 2 !(7Vi + A^ 2 )! 

X X! ^(^(X))— > X (t(Nx)) ^(^(JVi+X))— i^JVi+JV,))- (1-4) 

When {/j} is an orthonormal basis of fj, then {/^ V • • • V fy N }ii^.—^.i N is an orthogonal basis 
of fj^. Note that by definition 

/ v • • • v / = VM / ® • • • <g> /. 

AT times 

1.2. Creation and annihilation operators 

For every f E $), we define the creation operator a^(f) on := ^N^l^^/s C .F a / S by 

N r fiiV4 



requiring at{f)S)* /s C i}^ 1 for all TV ^ 0, with 



f Aip for fermions, 
V ^ for bosons. 
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By linearity, a'(f) can be defined as an operator on T^J S - Note that if {/iji^i is an or - 
thonormal basis of f), then {IlfcLi a Hfik)^}ii<---<iK, k^o is an orthonormal basis of T a and 
{]lfc=i at (/»fc) fi }»i<-<*jr, a:>o is an orthogonal basis of J" s . 

Similarly we define the annihilation operator a(f) by requiring a(f)Sj^ s C fi^fs 1 ^ or au 
TV Si 1, a(/)fi = and 

VV>eijf /s , (a(f)ip) (x 1 ,...,x N - 1 ) := VN I f(x)ip(x,x 1 ,...,x N - 1 )dx. 

jR d 

It can be verified that a(/) is the adjoint of a'(f) on J^J a - 

v*,*' e (*, t(/)*% a/< = (a(/)*,*')^ a/ .- 

In the fermionic case the creation and annihilation operators satisfy the so-called Canon- 
ical Anticommutation Relations (CAR): 

a{g)a\f)+a){f)a{g) = {g,f)tr a , 
ot(/)at( ff ) + ot( fl )ot(/)=0, (1.5) 
a(/)a(fl) + o(ff)a(/) = 0. 

These relations are satisfied on but it is deduced from the CAR that ||a^(/)|| = 
ll a (/)ll — I/I^i hence that a^(/) and a(f) can be extended to bounded operators on the 
whole fermionic Fock space F a . In the bosonic case, the creation and annihilation operators 
satisfy the so-called Canonical Commutation Relations (CCR): 

a(g)aHf)-aHf)a(g) = (g,f)tr 3 , 
at(/)at( 5 )-at( 5 )at(/)=0, (1.6) 
a(f)a(g) - a(g)a(f) = 0. 

These relations are satisfied on Tf n . Now a(f) and a'(f) are unbounded operators. However, 
they are bounded on Ff N (with values in J 7 ^ 1 ) for every fixed N. 



1.3. Observables 

We now define operators and quadratic forms on T a i s . The most important one is the 
so-called number operator which equals N on any Sj^ g : 

AA:=0A. 

This operator is unbounded on F a / S and its maximal domain is 

V(M):=U = W°,il> 1 ,...)er : J2 N2 ¥ N t- <0 °}- 
More generally, for every (densely defined) self- adjoint operator A on Sj, we may define 



by 



JV^I \i=l J 
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the operator on J- a / s . When A is bounded from below, the domain of X^iLi A Xi is simply 



f\" V(A) C io^y s ; in the general case, YU=i ^. i s essentially self-adjoint on f\" V(A) C 
^a/s> see E2- The operator A is self-adjoint on the domain 



N . 















2 


X>(A) := < 












< 00 > 

















In the literature, the second quantization A of A is often denoted by J^i Ai or by dT(A). 
Note that A/" is the second quantization of the identity on f). 

The operator A can be expressed in terms of creation and annihilation operators. Let 
{fi]i^i be an orthonormal basis of A, with /j £ T>(A) for every i ^ 1. Then we have (both 
in the fermionic and bosonic cases) 



A = J2 ^{Afj) a(fj) = J2 A v ot (/0 a (^), with A ij = (f h Af J ) i 



(1.7) 



The above series are well-defined when restricted to any T>(A) C C ^" a /s an d they 

coincide with which is the correct interpretation of the (formal) equality (II. 7[) , 

Applying this to the number operator, we obtain: 



(1.8) 



Similarly, we can associate to any two-body operator W : S) 2 a / S — > ft^u an operator W 
on Fock space, defined by 



W:=o ®o0 J2 W ^ 



N^2 \l^i<J^N 

where Wij denotes the operator W acting on the variables Xi and Xj but not on the other 
variables. We do not discuss problems of domains for shortness. As for one-body operators, 
the second quantization W in Fock space of a two-body operator W can be expressed in 
terms of creation and annihilation operators as follows: 



W: 



£ Wi i ,Ma t (/i)ot(/ i )o(//)o(/ fc ), 



with 



W, 



(fi A fj, Wfk A $1)^2 
(fiV fj,Wf k V f ( )^ 



(fermions) 
(bosons). 



(1 + 5ij)(l + 8kt) 

Note the normalization factor (1 + 5ij)(l + She) = \\fi V fj\\ 2 \\fk V fe || 2 for bosons. 
In particular, for an iV-body Hamiltonian of the form 

JV 



H V (N) :=£ 



-A, 



W{x k -x e ), 



(1.9) 



(1.10) 



(1.11) 
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with the convention that H v (1) = — A/2+V and H v (0) = 0, the corresponding Hamiltonian 
in Fock space defined by M v := H v (N) can be expressed as 

H V = ]T hij aHfi) a(Jj) + Yl Wi ^M <ftfi) "Hfj) «*(/*) <h) (1.12) 
with Wij^ki as in (|1.10[) and 

Remark 1.1. Physicists rather prefer to use the creation operator <p'(x) of a particle at 
formally related to the 'smeared' operator a'(f) by 

a f (/)=/ f(x)^(x)dx, ^{x) = YT(x)a\f i ) 

JRd i>l 

where is any orthonormal basis of L 2 (R d ). The formula (|1.12p can then be rewritten as 
follows: 

H= / ( ^V(^ f (a;) • V<p(z) + V(x)<p i (x)p(x) J cfe 

+ - / / T^(.T-y)^ t (x)^ t (y)^(y)^(x)dxdy. (1.13) 



1.4. States, density matrices 

A (mixed or normal) state on a (separable) Hilbert space X is a non- negative trace-class 
self-adjoint operator T 6 ©i(X) such that Tr(r) = 1. A pure state is an orthogonal projector: 
r = |^E')(5 , |. By the spectral theorem, any state is a convex combination of pure states: 

r = Y n i l*i)(*il> where n * > and Ui = L 

Even when the system is expected to be in a pure state, mixed states are very important 
tools that we use all the time. 

We always use the word 'state' for mixed state and only make comments related to a 
more general notion of states (a positive and normalized linear form on a C*-algebra 0E]). 
We denote by 

s{x) := {r = r* > o : TheCr) = 1} 

the convex set of all states on the Hilbert space X. The natural topology on S(X) is that 
induced by the strong topology of &i(3t). The set S(3t) is convex but it is not closed for 
the weak-* topology of @i(3£). Indeed we have in general that if T n — T weakly-* with 
{r„} C S(3C), then r = T* > but 

Tr x (T) ^liminfTr3e(r„) = 1 

71— >-00 

which is the operator version of Fatou's Lemma [58]. However it is known [121 [58] that if 
Trx(r) = 1 then the convergence is strong: |r„ — r|| ei ^ — > 0. The fact that a weak-* limit 
of a sequence of states is not always a state is a disease that will be repaired in Section [21 
when we introduce the geometric topology. 
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For a state T on the fermionic or bosonic Fock space J- a / s , we define the density matrix 

a/8 J a/s 



|T]0,9) .ff ^sf by the relation 



( 51 o ... o g p , ^M/, o • • • o = TV^ /s (V a f (/i) ■ • • a\f q ) a(g p ) ■ ■ • a( 5l )) (1.14) 

where o = A for fcrmions and o = V for bosons. When p = q we use the notation [T]& for 
the usual p-body density matrix of T. Note that [r)(°> = Trjr , (r) = 1 by definition. 

Remark 1.2. // T commutes with the number operator M, that is T = Q) n>0 G n with 
G n : ->■ , then it holds \T]M = forp^q. o 

Remark 1.3. We may define by the same formula the density matrices [T]( p ' q " > of any trace- 
class operator T (not necessarily self-adjoint and non-negative). o 

For fermions the creation and annihilation operators are bounded and (|1.14|) always 
properly define the operators [r]( p,<J ). For bosons, however, assumptions on V are needed to 
make (I1.14| meaningful. In the following we almost always consider states on the truncated 
Fock space F^/s ^ or wmcn (|1-14|) makes sense, as we explain below. 

Any state G on the TV-body space f)^/ s can also be seen as a state on the Fock spaces 



•^a/s an< ^ -^a/si by extending it to zero on sectors of k particles with k ^ N . A calculation 
shows that the kernel of [G]^ is given for p = 0, N by the well-known formula: 

[G]^(xi, ...,x p \x' x , ...,x' p ) 

= ( )[ dy p+1 --- [ dy N G{xx,...,x p ,y p+ \,...,y N ;x' 1 ,...,x p ,y p+ \,...,y N ). (1.15) 

\P / JR d JR d 

Saying differently, it is obtained (up to a constant) by taking a partial trace of G with respect 
to N — p variables. In particular it holds Tr^[G]^ = Q. If p > N + 1, then [G]& = 0. 
If r is any state on the truncated Fock space , then [r]( p ' ? ) = 0i£p^N + lor 

q ^ N + 1. Furthermore all the [r]^' 13 ) are trace-class operators, as stated in the following 
fundamental result. 

Lemma 1.1 (Density matrices of states on J 7 ^). (i) For all ^ p,q ^ N and all 

state r € S(J-^), the density matrix [r]' p ' 9 ' is trace-class: 




Furthermore the map 

r a C ( 1 

a/ 

is continuous. 



r e s [T]M e 6l (*l /s ,tf a/s ) (i.i7) 



(ii) States on are fully determined by their density matrices: ifTi,T2 € S^J 7 ^) are 

such that \Ti]M = [r 2 ](P<d for all0^p,q^ N, then T 1 =T 2 . 

The bound (II . 16[) is certainly not optimal and it is only provided as an illustration. 
It is a well known general fact that (regular) states are fully determined by their density 
matrices [B]. Our proof below is based on the explicit relation (11.20)) between the density 
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matrices [r]*^' 9 ) and the state T, when the latter is in S{F2/ S )■ These relations are useful 
in practice. 

Note that the linear map in (| 1 . 1 T|) is not weakly-* continuous. If for instance {ip n } is 
an orthonormal system of and T n — \^f n )(^f n \ with \I> n = ipi A ip n £ jrj|, then r„ — 
weakly-* in 6i(-7\^ s ) but [r,^ 1 ) — ^ |</?i)(</?i| weakly-* in 6i(.fj). Indeed, the purpose of the 
next section is precisely to introduce and study a weak topology that renders all the maps 
in (|1.17[) weakly continuous. 

We now provide the proof of Lemma 11.11 



Proof. We start by proving that for any state T, it holds [r]' p ' 9 ' £ ®i(^a/ s >^>a/s) ^ or 
all < p,q ^ N. We introduce the matrix elements G mn — IL„rn„ : #"/ s — > ?)™/ s 
where II n := t^ n j(Af) is the orthogonal projector onto Sj^/a- Since T £ (Si^f^), we have 
G mn £ &i(Sj2/ s ^T/s) for all m,n ^ N. 

It is easy to see from the definition of the density matrices that [G mn Y P ' q ' > = except 
when to — p = n — q ^ 0. A calculation shows that, in terms of kernels, 




(1.18) 

Since the partial trace of a trace-class operator is itself trace-class, we conclude that 
[G mn ]( p ' 9 ) is trace-class for all ^ p, q ^ N, and that 



[G 



||G mil | Si ^Ti 




where we have used that ||II m riI„| ei ^5 |r| 6l = 1. The continuity in the trace norm is an 
obvious consequence of the continuity of partial traces. 

Let ^ p, q ^ N and recall that only the matrix elements G mn such that to— p = n—q ^ 
contribute to For instance [T]( N >® = [G Nk ]( N ^ = G Nk and [T](MO = [G kN ]^' N ^> = 

GkN for all ^ k ^ N. Indeed the following holds for all ^ p, q ^ N: 

min(AT— p,iV— q) 

[T]M= J2 i G p+^+^ q) (1-19) 

which implies () 1 . 1 6[) . If we think of the density matrices [r]^' 9 ^ as being given, the previous 
equation (|1.19[) is a triangular system which allows to find all the G mn by induction. Inverting 
this system leads to the following formula: 

min(./V— m, N— n) , , 

G mn = [lf m ' n > + i- 1 ) 3 [r] (m+J ' n+j) . (1.20) 

i=i 



This shows that states on J^, a are uniquely determined by their density matrices. □ 



For to = n, (II .20)) may be written 

N-m , L 

G mm = [T]^ + ( I ) Tr m +i^„ l+J [T]^ m +^ (1.21) 



i=i 



m 
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where Tr m+ i_j. m+J - denotes the partial trace with respect to the j last variables. For instance, 
we have [T]( N > = Gnn, and 

G N - 1N - 1 = [r)( N -V-NTr N [r]W 

which follows from the fact that 



Remark 1.4. Lemma \ 1 . 1\ is not true as such on the set iS(J r a / s ) of states on the whole Fock 
space. In general, we have ^ and 

which is finite only under appropriate assumptions on T. The off-diagonal density matrices 
[p](p.?) are i n general only Hilbert- Schmidt when all the |T]( p ' are trace-class. o 

Remark 1.5. We say that a family of operators {T" l }^ =0 with T m G &i(S) m ) is J 7 ^- 
representable when there exists T £ with [T,Af] = such that r( m ' m ) = T m for 

all m — 0,...,N. Using Formula (jl.21jl . we see that -representability is equivalent to 
having T° = 1 and 

N , .* 

Vm = 0, N, T m + {-l) 3+m (I) Tr m+1 ^- > 0. 

j— m+l 

The case of states which do not commute with Af is more involved. o 

In this section we have introduced Fock spaces and creation/annihilation operators for 
indistinguishable fermions or bosons. When working in the truncated space defined in 
(|1.2[) . the statistics of the particles does not make a big difference. To simplify notation, we 
now write $) p , J 7 , J 7 ^ 1 ^ , etc, without specifying the considered statistics, except for results 
which are specific to bosons or fermions. 



2. Geometric convergence 
2.1. Definition and properties 
2.1.1. Definition 

We define a weak topology on states in ^(J 7 ^), induced by the weak-* topologies of all 
the density matrices [r]^' 9 -*: 

Definition 2.1 (Geometric topology &: convergence). We define the geometric topol- 
ogy T on S(J 7 '^ N ) as the coarsest topology such that the maps 

reS{^ N ) h-> U, [r]^y\ (2.1) 

\ /ftp 

remain continuous for all ip') € f) p x fy q and all ^ p, q ^ N . 

Let {r„} be a sequence of states on , and T be a state on J 75 *^. The sequence {F„} 
is said to converge geometrically to F if 

lim U, [r„]^y\ = U, [T]^^) (2.2) 
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for all (ip, tp') £ S) v x S) q and all ^ p,q ^ N. We use the notation T n T. 

g 

Note that, when it exists, the geometric limit V is uniquely defined since T <= S(J-^ N ) is 
characterized by its density matrices [r] , by Lemma 11.11 

We give several examples right after the following result which is an immediate conse- 
quence of Lemma 11.11 



Lemma 2.1 (Elementary properties of geometric convergence). 

(1) The geometric topology T is coarser than the usual norm topology. If L„ — » T strongly 
in (3i(J r ^ A ') 7 then T n ^ g T geometrically. 

(2) We have T n -^ g T in , if and only if \T n ]M [T]^ weakly-* in & 1} for all 
0^p,q^N. 

Proof. The first assertion follows from the (strong) continuity of the maps T H> [r]( p ' 9 ) 
for all ^ p, q ^ N, as stated in Lemma 11.11 The second assertion is a consequence of 
the uniform trace-class bound (|1.16[) on all the density matrices and of the Banach-Alaoglu 
Theorem in 6i(55 9 ,^ p ). □ 



Let us emphasize that the geometric limit L of a sequence of states is, by definition, always 
a state, that is it must satisfy Trjr(L) = 1. Contrarily to the usual weak-* convergence on 
6i(J r ^ Ar ), there is never any loss in the trace-norm when L„ — ^ V. If in the geometric limit 
some particles are lost, then L lives on spaces with less particles in J 7 ^. If all the particles 
are lost, then we have L = |0)(f2|, the vacuum state in J 7 ^. 

We now provide examples of sequences {L„} which geometrically converge but do not 
strongly converge to a limit T. Our claims can be verified by computing the density matrices 
[r„]( p ' 9 ) and checking their weak-* convergence towards [L]( p ' 9 ). 

Example 2.1. Let {(p n } be an orthonormal basis of £). Define a sequence of two-body 
fermionic wavefunctions by := p\ A p n , with associated state in T^ 2 denoted by L„ = 
© © |*„)(*„|. It holds T„ weakly-* and r„ -^ g © |<£i)(^i| © geometrically in 
J 7 ^ 2 . The geometric limit T describes a system composed of only one particle, in the state 
<pi. The other particle in the state <p n has vanished in the limit. o 

Example 2.2. Even when L„ is a pure state for all n, the geometric limit L is not always 
a pure state. For instance if ^„ := p\ A p 2 n with p x n = cosapi + sin a ip n and p 2 n = 
cos /3 f2 + sin (3 <p n +i, then the corresponding state F n = OffiOffi |\E , „)(\I'„| on J^ 2 converges 
geometrically to 

r„ L = (sin 2 a sin 2 /?^ (cos 2 a sin 2 /3 | | + cos 2 /3 sin 2 a \ip2){v>2\) 

© (cos 2 a cos 2 /3 \ipi A tp2){<-P\ A ■ 

On the other hand, we have 

r„ 8 0® cos 2 a cos 2 /3 |</?i A </? 2 )(</?i A tp 2 \ 

weakly-* in S^J 7 ^ 2 ). o 

Example 2.3 (Hartree states). For bosons, a Hartree state takes the form ^ = (p® ■ ■ ■ (g) 
P € io^ where = 1- Assume that {<p n } is a sequence of normalized functions in .f), with 
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ip n ip weakly. Let T n = © • • • © \(<Pn)® N ){(.<Pn)® N \ be the associated iV-body state in 
S(Tf N ). Then it holds 



Example 2.4 (Coherent states). For bosons, coherent states are defined by the formula 
Tf ;= W(f )\fl) £ F s where W(f) = exp(a i '(/) — a(/)) is the Weyl unitary operator (/ is any 
vector of the one-body space S) 1 ). The latter satisfies the following interwinning relations 

W(fy(a(g)-(g,f))w(f) = a(g), W(f)* (a\g) - (f,g))w(f) = a\g). (2.3) 

The density matrix [Y of a coherent state T f = W(f) W(f)* is [T f ]^ = 
I /® p ) (/® 9 1 . Consider a sequence {r * } of coherent states with {/„ } bounded in Sj 1 , such that 
/„ / weakly in ft 1 . Then T fn ^ g T f geometrically in the sense that [T fn ]^ ^„ [Tf]^ 
weakly-*, for all p, q ^ 0. Note that coherent states do not live on any truncated Fock space 
J 7 ^ 1 ^ , hence Definition 12.11 has to be generalized in an obvious fashion on the whole Fock 
space J ' . o 

Example 2.5 (Hartree-Fock(-Bogoliubov) states). For fermions, there is a subclass of 
states which are fully characterized by their one-body density matrix [T]' 1 ) and their pairing 
density matrix [r]^' - 1 (if they commute with TV, they are only characterized by [FJ^ 1 )). 
These states are called generalized Hartree-Fock states [4] or Hartree-Fock-Bogoliubov states 
(when [r^ 2 ' - 1 ^ 0). Here 'fully characterized' means that any density matrix [r]( p ' 9 ) is an 
explicit function of [r]^ 1 ' 1 ) and r^ 2 '°\ given by Wick's formula, see Eq. (2a. 11) in [4]. When 
Tr^ ([TV,]' 1 )) is uniformly bounded, it is easily seen that geometric convergence of generalized 
Hartree-Fock states is equivalent to the weak-* convergence of [r n ]W and of [F n ]( 2 ' ). The 
geometric limit is always a generalized Hartree-Fock state. 

Note that if [r]W has an infinite rank (but a finite trace), then the corresponding Hartree- 
Fock state r does not live on any truncated Fock space J r ^ fe . However, geometric convergence 
can be understood in the same fashion as in the previous example. o 

Example 2.6. Let r be any state on and let U{t) = e~ itT (with T = - A/2) be the 

unitary free evolution on the one-body space $y, of a non-relativistic particle. Let 



be the unitary evolution of the second quantization of T on the truncated Fock space T 



The state T(t) := V(t)ToV(t)* is the unique weak solution to the Schrddinger-von Neumann 
equation 




It is clear that the convergence is strong if and only if ||y|L = 1. 



O 






18 



Mathieu LEWIN 



Then 

r(t) |fi)(fi| asi->±oo. 

Indeed, we have 

VO < p,q ^ N, [T{t)]^ q) = U(t) <g> • • • <g> I7(i) [r ] (p ' 9) U(t)* ® ■ ■■ ® £7(t)* 



P 9 

which tends to weakly-* in ©i, when (p, q) 7^ (0, 0). The same holds if U(t) is any unitary 
family satisfying U(t) weakly as t — > ±00. o 

After these examples, we now make some fundamental remarks about the notion of 
geometric convergence. 

Remark 2.1 (Geometric convergence is a C*-algebra concept). The geometric topol- 
ogy is the restriction to of a well-known weak topology arising in C* -algebra theory, a 
fact that we will need in the proof of Lemma \2.2\ below. 

For fermions, an equivalent way of formulating (|2.2I) is, by the definition (| 1 . 14[) of density 
matrices, 

VA G A, lim Tr^(r„A) = Tr^(IM) (2.4) 

n— y 00 

where A is the C* -algebra generated by all the (v(f) with f any vector infi. Therefore 

for fermions the topology T is nothing but the usual weak-* topology of states on the CAR 
algebra A, restricted to states of the truncated Fock space F^/s ■ For bosons, the same holds 
true with A being the CCR algebra, generated by the Weyl operators of Example \2.4\ 
Note that we have T n T for the weak-* topology of ©i(J r ) if and only if 

VK G lim Tr^(T n K) = Ti^TK) 

n— too 

where we recall that K.(J-) is the algebra of compact operators. In both the fermionic and 
bosonic cases, the CAR/CCR algebra A does not contain any nontrivial compact operator: 
A(~\IC(J-) = {0}. Geometric convergence is thus a priori not related to the usual weak-* 
convergence and it is possible to have F n — ^ T with Trjr(T) = 1 whereas T n — weakly-* 
in &i(jF), like in the previous examples. o 

Remark 2.2. If T n commutes with the number operator M for all n, \T n ,N] = 0, then 
[r„p>«) ee for allp^ q and it is easy to verify that the geometric limit T of {T n } must 
also commutes with M . o 

Remark 2.3. A similar definition of the geometric topology and convergence can be pro- 
vided if the system contains several species of particles. One introduces the density matrices 
ppj(pi,...,P«,,gii...i5/) w fi ere p i an d g i respectively count the number of annihilation and cre- 
ation operators of the species i (bosons or fermions). One works in the truncated Fock space 
jr^N 1 ,...,N k correS p 0n( H n g t having at most Ni particles of species i. o 

2.1.2. Compactness results 

The following result is very useful in practice. It allows us to work with weak limits of density 
matrices while being sure, at the same time, that the limits arise from a state T. 

Lemma 2.2 (Geometric compactness of S(J-^ N )). The set of states S(J-^ N ) on J 7 ^ 1 ^ 
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is (sequentially) compact for the geometric topology T: every sequence of states {r„} C 

S{F^ N ) has a subsequence which converges geometrically, T nk — T. 

s 

Proof. This result immediately follows from well-known facts in the theory of C*-algebras 
(recall Remark 12. By the Banach-Alaoglu Theorem, any sequence of states {r„} c 
iS^^) on the CAR (resp. CCR) algebra A generated by the creation operators (resp. 
Weyl operators), has a weakly-convergent subsequence in the sense that for every A G A, 
one has Tr(r„ fe ^4) — > uj(A) where w is a positive normalized linear form on A, [5]. Since 
r„ lives on the truncated Fock space J-^ N for every n, it has a uniformly bounded average 
particle number, hence its weak limit uj must be a normal state [6]: there is a T € ^(J 7 ) such 
that lo{A) = Tvjr{YA) for all A. Since [r]( Ar + 1 - Ar + 1 ) = 0, it is easy to verify that V must also 
live on and the result follows. □ 

Remark 2.4. Up to extraction of subsequences, one can always assume that [r n ]^ p ' 9 ' 1 
T( p ' 9 ) weakly -* in &i(S) q ,f) p ). The matrix elements G m ,ri of the limit state T are then 
uniquely determined from the operators by Formula (|1.20p . What is more subtle is the 

fact that the so- obtained T is really a state, that is T = T* ^ 0. o 

Remark 2.5. Lemma \2.2\ can obviously be extended to sequences of states {r„} on the 
whole Fock space T which satisfy a uniform bound of the form Trjr(A/T„) ^ C . o 

The following result says that the total number of particles in the system cannot increase 
under geometric convergence, and that there is strong convergence if and only if no particle 
has been lost. 

Lemma 2.3 (Average particle number and strong convergence). Let {T n } be a 
sequence of states in S(J-^ N ) and T € S(J-^ N ) be a state such that T n —^ g T . The average 
particle number is lower semi- continuous: 

Trjr(A/T) sC liminf Trjr (A/T„ ) . (2.5) 

n— > oo 

Furthermore, i/ lim„_>. 00 Trjr(J\fT n ) — Trjr(J\fT), then T n T strongly in &i(JF^ N ). 

Proof. Let us recall that Tr^(A/T) = Tr^rp), hence, since [r n ]W [T]W weakly-* in 
©i(f)) by Lemma [2~Tl we have 

Tr^A/T) = TiY,[r] (1) sC liminf Try, [r„] (1) = liminf Tr^fA/TJ. 

Another proof consists in writing that Trjr (r n YliLi a Hfi) a (/i)) ^ Tr^r (T n Af) by (|1.8p . It 
then suffices to pass to the limit first as n — > oo and then as K — > oo. 

The proof that conservation of the average particle number implies strong convergence 
requires a bit more work. We start with a sequence of A-body states, that is T„ = Offi • • - ®G r ' 

commutes with J\f: 



where G™ € S(Sj N ). We assume that T n T in . From RemarkES we know that T 



/Goo 
\ Gnn/ 
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The assumption that 

N 

N = lim Tr^(r„) = Tr^(WT) = V k Tr^ (G kk ) 

n— ^oo L — * 

k=Q 

together with the fact that ^2 k=0 T r ^ fe (Gfcfe) = 1 since G is a state, imply that Gkk = for 
all k — 0, TV— 1 and Tt^n (Gnn) = 1. However, we know that Gnn is the weak-* limit of 
G™ in &x(f) N ). Therefore Tt^n(G nn ) = 1 implies that G" ->■ G NN strongly in 6i, by the 
reciprocal of Fatou's Lemma for trace-class operators (see [HI [S3]), and the result follows. 

We now come back to the general case. Let F„ ^ g F be an arbitrary sequence which 
converges geometrically in J r ^ N , such that Trjr(A^F) = lim„_ i . 00 Tr^A/T). We denote by 
G ke the matrix elements of F n and introduce the auxiliary state 




obtained by retaining only the diagonal of T n . It is easy to check that f n — * a T, the diagonal 
of r. We first prove that F„ — > F strongly. Indeed we may write 

N 

f n = ^ tk&kk 
k=0 

where t k = Tr^k (G kk ) and G kk = G kk /t k is a state on S) k (with an obvious convention when 
q = 0). We have G n kk ^ g G kk for all k = 0, N and f = J2 k =o with t k = lim^^ q 

(up to subsequences). Our assumption means that 

N N 
k=0 k=0 

However by (|2.5|) . it holds Trjr(J\fG kk ) ^ k for all fc, hence the previous equation means 
that Trjr(NG kk ) = k for all k = 0, N such that t k ^ 0. As we have shown in the previous 
paragraph, this implies that G kk — > G kk strongly in ©i(f) fc ). When t k — 0, we have simply 
G2 k — > strongly. This eventually shows that r„ — > F strongly. 

We now conclude that r„ — > T strongly. Indeed, we have F„ — 1 L' weakly-* in &i(F^ N ) 
and we know that the diagonal of F„ converges strongly, hence in particular Tr(F') = 1. By 
the reciprocal of Fatou's Lemma [12j [58] , this implies that r„ — > T strongly, which ends the 
proof of Lemma 12.31 □ 



2.2. Application: HVZ theorem in the lower semi- continuous case 



In this section, we illustrate the use of geometric convergence on the very simple example 
of a many-body system with a non-negative two-body interaction. Our example covers the 
celebrated case of atoms and molecules. 

We consider the following many-body Hamiltonian 



H V (N) :=£ 
j'=i 



V( Xj 



W(x k - x t ) 



l<k<e<N 



(2.6) 
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on L^ s ((R d ) N ). Since in practice W is fixed (it is a characteristics of the studied particles) 
whereas V is an external field that can be varied, we only emphasize V in the notation of 
H V (N). We choose any statistics (bosons or fermions) for our particles. The spectrum of 
H v (N) depends on this statistics but our results are stated the same in both cases. 

We assume that W is even and that the two real functions V and W can both be written 
in the form 53i=i fi with fi G L Pi (R d ) where max(l, d/2) <pi < oo or pi = oo but — > at 
infinity. These conditions ensure that (1- A) _1 / 2 V(1- A) -1 / 2 and (1- A)-^ 2 W(1~ A)" 1 / 2 
are compact operators. Then, by the KLMN Theorem [5T|, H V (N) has a unique self-adjoint 
realization in the TV-body space L 2 ^ s ((M. d ) N ) with quadratic form domain H*,((Mr) N ). 
More precisely, for every < e < 1, there exists a constant C = C(N, e) ^ such that 




(1-e) > -A Xj -C^H v (N)^(l+e) > -A Xj +C (2.7) 



in the sense of quadratic forms on L^ s ((M. ) ). In this section we will make the assumption 
that the interaction is repulsive, that is 



W > 0. 



The general case is treated later in Section 13.21 

Example 2.7 (Atoms and molecules). For atoms and molecules in which the electrons 
are treated as quantum particles whereas the nuclei are considered as fixed pointwise classical 
particles (Born-Oppenheimer approximation), we have in atomic units, on L 2 ((M. 3 ) N ), 

M 

where R m and z m are the positions and charges of the nuclei. The functions V and W 
are respectively the Coulomb attraction potential induced by the nuclei, and the Coulomb 
repulsion between the electrons. o 

The second-quantization of H v (N) is the Fock Hamiltonian 

M v :=Oe0ff y (fe) 

which we restrict to the truncated Fock space J 7 ^. The energy of the system in the state 
r G S (J 7 ^) reads, using (|1.12l) and the definition (II . 14)) of the one- and two-body density 
matrices [r]W and [r]( 2 ^: 

£ V (T) := Trjr(m v T) 

= Tr L2(Rd) (Y-Ia + f) [r]«) +Tv Ll/siWLdxMd) (V[r] (2) ) • (2.8) 

By (|2.7[) . the energy is well-defined for states T G S(J-^ N ) such that 
Tr^ ( Ji/W 2 ) = Tr fl ( (-A^r^-A) 1 / 2 ' 



Tr fl2 ({~A X + -Aj / ) 1 / 2 r' 2 ^(-A a , - Ay) 1 / 2 ] < oo 



N-l n <-i* 

When the previous kinetic energy term is infinite, we can let £ v (T) := +oo. 
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One difficulty of many-body systems is the lack of weak lower semi-continuity (wise) 
of the quantum energy ^ £ Sj H> (ty, H v (N)^. This was for instance pointed out by 
Friesecke, see Lemma 1.2 (iii) in [T!5]. Indeed if we denote by 



E V {N) := inf a (H v (N)) , E y (7V) := inf a css (H v (AT)) 



(2.9) 



respectively the ground state energy and the bottom of the essential spectrum, we usually 
have that T, V (N) < 0. This implies that for a singular Weyl sequence >]/„ — ^ it holds 
H v (N)^ n ) Y, V (N) < 0, showing that the energy is not wise. 
We now prove that, on the contrary, when W ^ the energy is lower semi-continuous 
for the geometric convergence which we have introduced in the previous section. 

Lemma 2.4 (Lower semi-continuity of the energy under geometric convergence). 

Assume that W ^ and let {T n } be a sequence of states in which converges geomet- 

rically to r. Then 



£ V (T) ^ liminf£ y (r n ). 

n— 7-oo 



Proof. Under our assumptions on V and W, it is easily verified that £ v is lower semi- 
continuous for the strong topology of &i(J-^ N ). We have to prove that the same holds for 
the geometric topology. 

When the kinetic energy of {T n } is not bounded, there is nothing to show by (|2.7|) . hence 
we may as well assume that 



Tr^ (V /2 T„ T 1/2 ) < C 



for a constant C independent of n (this is actually equivalent to assuming that each p-body 
density matrix has a bounded kinetic energy). Since we have by assumption |T ra ]( p ) — 1 [r]( p ) 
weakly-* in 6i, we deduce that the geometric limit T has a finite kinetic energy, hence a 
finite total energy. We now remark that 

£ v {T n ) = ^Ti L 2 {Rd} ((-A)[r„]W) + J ^ Vpr n +T rLl/ARdxMd) (w[T n }^ . 

where pr n {x) — [r n ]W(:r, x) is the density of the system. It is then a classical fact that 
Tr L , (Rd) ((-A)[r]«) ^ liminf Tr L2(Rti) ((-A)[r„]«) 



Tr 



T 2 



ldxRd) (W[T]W) ^ liminf Tr Ll/B(Vxm (w[T n ]W 



Vpmw = um / Vp\r„]W- (2.10) 

The first two claims follow from Fatou's Lemma for trace-class operators [57] (using W ^ 0). 
The last claim (12.101) is shown as follows. First the Hoffmann-Ostenhof inequality [24] 



/ |vVpF| 2 < Tv mRd) ((-A)[r]«) , 



(2.11) 
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implies that Jpr n is bounded in H 1 (R d ), hence we may as well assume that JpT„ 



weakly in H 1 ^) and strongly in Lf oc (R d ). Recall that V = J^f=i v j with v j € £ Pj (M d ) 
where max(d/2, 1) < pj < oo or Vj E L 00 (R rf ) and Vj — > at infinity. For d > 3, we write 

y?' ( x ) Pr„ (a;) dx ^ l^jfli« , i(M d \B(o 1 Ji)) l^r n |li"i(ii<i) ^ ^WnL v i{m d \B{<3,R)) (2-12) 

\x\^R 

where 1/pj + 1/qj = 1, hence 1 < qj < d/(d — 2). In the last inequality we have used the 
Sobolev injection theorem as well as the fact that y/pr n is bounded in H 1 (M. d ). On the other 
hand, by Rellich's theorem, we have a compact injection H 1 (B(0, R)) °-> L q (B(0, R)) for all 
2 ^ q < 2d/ (d — 2) which implies that 



lim / Vj(x) pr (x) dx = / Vj(x) pr{x) dx. 

n ^°°J\x\^R J\x\^R 

Together with f)2 . 12[) . this proves (|2 . 10|) . The proof is the same in dimensions 1 and 2. □ 

The following is a famous result for many-body systems: 

Theorem 2.1 (HVZ in the lower semi-continuous case). Assume W ^ 0. Then it 
holds E°(N) = for allN^O and 



E y (iV) = E V (N- 1). (2.13) 



In particular, E v (N) is an isolated eigenvalue if and only if 

E V (N) < E v {N-l) = mxa{E v (N -k)+E°(k), k = 1,...,N}. 

Remark 2.6. A similar result holds true if the system contains several kinds of particles 
(with possibly different interaction potentials) , with or without internal degrees of freedom. o 

Theorem 12.11 is due to Zhislin [66] , Van Winter [62] and Hunziker [25] . Simpler proofs 
were provided later when the so-called geometric methods were developed [16l [57] [54j [9] . 
The interpretation of (|2.13[) is that in order to reach the bottom of the essential spectrum, 
we have to provide a sufficiently large amount of energy to the system in order to extract 
at least one particle. The case of a general interaction W is treated later in Section |3~21 

Theorem l2.1l is essential when proving existence of ground and excited states. The bottom 
of the spectrum E V (N) is an isolated eigenvalue if and only if the HVZ inequality E v (N) < 
E v (N — 1) holds. Such an inequality can be proved by induction on N: admitting that 
E V (N — 1) < E V (N — 2), there is a ground state for E V (N — 1) and one can use this state 
to construct an iV-body test state to prove that E v (N) < E v (N - 1). 

For atoms and molecules (Example 12 .7j) . Zhislin and Sigalov [66] [67] have shown that 
there is a ground state as well as infinitely many excited states as soon as N < Z + 1 where 
Z = J2m=i z m is the total nuclear charge. The idea is that, with N — 1 electrons bound to 
the nuclei, any additional electron escaping to infinity sees a Coulomb interaction induced 
by a total charge Z — (N — 1) > 0. This potential is attractive at large distances and the 
desired inequality E v (N) < E v (N - 1) follows. 

We now turn to the proof of Theorem 12.11 

Proof. The bound (N) ^ E V (N — 1) is shown by building a convenient singular Weyl 
sequence, using a Weyl sequence for E v (N — 1). We do not elaborate more on this classical 
fact and we only explain the proof of the more complicated inequality T, V (N) ^ E V (N — 1). 
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First we note that since E v (N) < T, V (N) < ^^(iV - 1), the map N ^ E^AT) j s 
non-increasing. When V = 0, £7° (AT) ^ since W ^ 0, hence E n (N) = for all AT. 

Let now {^ n } C i5 w be a singular Weyl sequence for T, v (AT), that is such that (H v (AT) — 
T, v (N))^ n -» 0, ||*„|| = 1 and weakly in L 2 ((R d ) JV ). The corresponding pure state 

onJ^ isT„ := 0©- • -0©|*„)(* n | and it has a bounded energy, lim n ^oo £ v (T n ) = E V (N), 
hence a bounded kinetic energy by (|2.7p . Extracting a subsequence if necessary, we may 

assume by Lemma \2. 2 1 that F„ — 1 F geometrically. We write as usual 

g 

/G 00 

r = 

\ Gnn/ 

Recall that Gnn = [Gnn]^ N ' is the weak-* limit of |^ r „)(4 r n |, hence Gnn = since — 1 0. 
By Lemma \2. 41 we have 

N-l 

Z V (N) = lim E v {T n ) > E V {T) = V Tr^- {H v {j)G n ) 

n— foo *■ — * 

^ E V (j) Tr^ (Gjj) >E V (N- 1). (2.14) 

3=0 

In the second line we have used that Gjj ^ and that Sjlo 1 Tr^j {Gjj) — 1 since F is a 
state. □ 



Remark 2.7. Let {^ n } be a singular Weyl sequence for the bottom X 17 (N) of the essential 
spectrum of H v (AT) , like in the proof of Theorem \2.1\ Then, if E v (N — 1) < Y, v (N — 1) = 
E v (N —2) , it can be seen from (|2.14j) that its geometric limitT is a ground state of H v {N — 1) 
inLl /s ((M. d ) N - X )- ' o 



3. Geometric localization 

Localization is a fundamental concept of many-body quantum mechanics. In the seminal 
works of the end of the seventies [TH1 El [S71 SSI [Ml I2S] > the expression 'geometric methods' 
was used to denote the use of appropriate partitions of unity in configuration space. In 
this section we explain how one can lift a localization in the one-body space Sj 1 to the 
truncated Fock space J 7 ^ , following Derezihski and Gerard [13], and we relate this tool to 
the geometric topology defined in the previous section. 



3.1. Definition and properties 
3.1.1. Definition 

Here we explain how to localize a state T £ S(J-^ ). As already suggested in the introduc- 
tion, the localization of a pure one-body state ip € L 2 (M. d ) in a domain Del' 1 should be 
described by the state 

r x = (i-^Jx^| 2 ) ®\x<p)(x<p\ (3.1) 
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where x — 1 D ■ Note that the previous formula actually defines a state for every normalized 
ip G L 2 (R d ) and every function x such that ^ \x\ 2 ^ 1. This discussion suggests the 
following definition of localized states. 

Proposition 3.1 (Definition of localized states). Let B G B($j) be a bounded operator 
on Sj, such that ^ BB* ^ 1, and T G S(J-^ N ) be any state on J-^ N . Then there exists a 
unique state Tb G S(J-^ n ) such that 

[T B ] M = B ® ■ ■ ■ ® B \T]M fl* ® ■ ■ ■ ® B* (3.2) 

for all p, q ^ N . The state Tb is called the B-localization of T. 

Note that in general the localized state Tb is not a pure state, even when T is itself a 
pure state. The concept of localization of states in Fock space was first introduced for bosons 
by Derezihski and Gerard in |13j and generalized to fermions by Ammari in pQ. It is now 
a classical tool in Quantum Field Theory. It was recently used by Hainzl, Solovej and the 
author of the present paper, to prove the existence of the thermodynamic limit for quantum 
Coulomb systems in the grand canonical picture, see Appendix A.l in 23 . In this latter 
work, the strong subadditivity of the quantum entropy was also formulated using geometric 
localization. Although expressed in different terms, the definition of Tb in Proposition 13. II 
coincides with that of all these previous works. 

We now turn to the proof of Proposition 13.11 

Proof. A state satisfying (|3.2I) was constructed in [13l [IJ [23] , using the partial isometry 
/ ei)4 Bf® VI - BB*f G Sj © Sj and the fact that J"(#i © £2) ~ F(S)i) © The 
state is obtained by means of a partial trace with respect to the second Hilbert space. 
Uniqueness then follows from Lemma ll.ll □ 

Remark 3.1. The matrix components {Gmn}m n=0 of the operator Tb can be expressed 
using Equation (|1.20j) as follows 

min(N—m,N~n) , , 

qB _ ^®m|-pj(m,n)^*^®n _|_ (~ l) 3 _B®( m +-?) [F] ( m+ 3> n+ 3) (£}* 

3=1 

(3.3) 

The verification that the so-obtained operator is a state (T b = (Tb)* OJ uses the 
CCR/CAR algebra A in a similar way as in the proof of Lemma \2.2l o 

Remark 3.2. If B 1 and B 2 are such that < B k B* k < 1, then (B 2 B 1 )(B 2 B 1 )* = 
B 2 B\B* X B* 2 ^ B2B2 ^ 1. It is then clear from the definition that (Tb 1 )b3 = TsaSi- 

We now illustrate Propostion 13. II by several examples of localized states. 

Example 3.1. We have for all state Ti = T and To = |0)(fi| (the vacuum state), corre- 
sponding to having, respectively, B — 1 and B = 0. If Lp E f) 1 and T — © \ip)((p\, then 
T B = (1 - ll-B^f) © \B<p){B(p\, as in (£]}. o 

Example 3.2. If U is a unitary operator on £j\ then (T) a = (1 ©[/©••• © f/®^) T (1 © 
C/*ffi---ffi([/*)^ Ar ). o 
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Example 3.3 (Localization of iV-body states). Let G E S(S) N ) be an iV-body state 
and r = © • • • © G € ^(J 7 ^^). A simple calculation based on (I3.3[) shows that 

T B = Gf © • • • © G% 

where 



= (J^j Tr k+1 ^ N (V fe © y/l - BB* MN k) G (B* ) m ® \fl - B B '" ) C-U-i 



■®(N-k) 



with Trfc_|_i_j.Ar denoting the partial trace with respect to the N — k + 1 last variables. More 
explicitely, if G = |*)(*| and s$ xi x ) < 1, then 

G^(i 1) ...,^ari > ...,4)=(^)lI^iM^) /•••/ ft (1-X 2 fe))x 



x ^{xi, ...,x kl z k+1 , ...,z N )^(x' 1 , ...,x' k ,z k+1 , ...,z N )dz k+1 ■ ■ -dz N . (3.5) 
We see from (ET41 that it holds 



TV (Gf ) = Tr^_* (G£T k BBm ) 



(3.6) 



The relation (|3 .6[) will play a very important role later and it may be considered as one of the 
basic tools of the geometric methods for many-body systems. For B = trj(x), it essentially 
means that the 'weight' in the fc-particle sector of the localized state in a domain D is equal 
to that in the (N — fc)-particle sector outside D. o 

Example 3.4 (Hartree states). Let V = 0© • • • © \ip® N ){ip® N \ e Ff N be a Hartree state 
as in Example 12.31 Then 



k=0 



Example 3.5 (Coherent and Hartree-Fock-Bogoliubov states). If Tf is a coherent 
state like in Example 1 2. 4[ then (r j)b = Tb/- If T is a Hartree-Fock-Bogoliubov state like in 
Example 12.51 with one-body density matrix [T]W and pairing density matrix [r]( 2 '°), then 
T B is the unique Hartree-Fock-Bogoliubov state having B[T]^B* and (J5 © B) [r]( 2 -°) as 
one-body and pairing density matrices. In Example 14.21 below we detail the case of pure 
Hartree-Fock states. o 



3.1.2. Convergence results 

Let us now turn to some useful applications of geometric localization. We start by showing 
that the localization map T i— > Tb is continuous with respect to the geometric topology. 

Lemma 3.1 (Continuity of geometric localization). Let {T n } be a sequence of states 
in S(F^ N ) which converges geometrically to a state T € S(J-^ N ), T n T. Let B G 
yB(f) 1 ) be such that ^5 BB* 1. Then the associated sequence of localized states converges 
geometrically: (T u )b —^g Tb- 

Similarly, if B n is a sequence satisfying $5 B n (B n )* ^ 1, B n — > B and (B n )* — > B* 
strongly (that is B n x — > Bx and B* n x — > B*x strongly in fj 1 for any fixed x € then it 
holds (T n ) Bn T B - 
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Proof. When T n ^ g T, that is \T n ]M -S, [T]M for all p, q ^ N, we have that 
[(r„)B] (p ' 9) = B®P[T n ]M(B*)^ converges weakly-* to B®p[T]^ (B*)®* . This is by defi- 
nition [r B ] (p ' 9) , hence it holds {T n ) B -^ g T B . The argument is the same when B n — > B and 
(B n )* -> B* strongly. □ 

The next lemma explains how localization can be used to convert geometric convergence 
into strong convergence. 

Lemma 3.2 (Local compactness). Let T be a non-negative self-adjoint operator 
on Sj 1 , and B be a bounded operator such that ^ BB* ^ 1. We assume that B is T 1 / 2 - 
compact, that is B(l +T 1 / 2 ) -1 € /C(f) 1 ). Let {T n } be a sequence of states in S(J-^ N ) which 

converges geometrically to a state T G S(J-^ N ), T n T . If 

g 

Tr^rV^igWT 1 / 2 ) <C 
for a constant independent of n, then (r n )s — > Tb strongly in ©i^J 73 *^). 

Proof. We have (T n ) B — geometrically by Lemma [3.11 and it remains to prove that 
the convergence is strong. It holds 

[(r„) B ] (1) = B[T n ]^B* = k(i + r 1 / 2 )[r n ]( 1 )(i + t^ 2 )k* 

where K = ^(l+T 1 / 2 )" 1 is compact by assumption. The sequence (l+T 1 / 2 )^] W (1+T 1 / 2 ) 
is bounded in &i(Sj), hence we have that (I + T 1 l 2 )[T] {1 \l +T 1 / 2 ) e & 1 (S) 1 ) and 

(l+T 1 / 2 )[r„]( 1 )(l+T 1 / 2 ) _^ (l+^/^T^Hl + T 1 ' 2 ) 

weakly-* in &\. It is well known that if A n — A weakly-* in 61(f) 1 ) and K is com- 
pact, then KA n K* — » KAK* strongly in 61 (fj). We deduce from the above calculation 
that [(r„)s] (1) -> [r B ] (1) strongly in ©i(fj). By Lemma O this shows that (L„) B -)• T B 
strongly. □ 

Example 3.6. If T n —*- g V geometrically in and the kinetic energy Tr ((— A)[r n ]( 1 ^ ) ) 

is uniformly bounded, then (T n ) x — > T x strongly in 61, for every localization function x( x ) 
of compact support (even tending to zero at infinity), since x(^)(l + I ^ *V|) _1 is always a 
compact operator. This can be viewed as a generalization to states in J 7 ^ 1 ^ of Rellich's local 
compactness in Sobolev spaces [33]. o 

The following is simple consequence of the previous result with T = 1 . 

Corollary 3.1 (Compact localization). Let {T n } be a sequence of states in iS(J r ^ Ar ) 
which converges geometrically to a state T £ S(F** N ), L„ — x g T. Then (T n )K (T)k 
strongly in ©i(J r ^ Ar ) for every fixed compact operator K such that KK* 1. 

Localization may also be used to approximate a given state by simpler states (for instance 
finite rank states, see Section |4}. 

Lemma 3.3 (Approximation by localized states). Let {B n } be a sequence of bounded 
operators in f), such that ^ B n B^ ^ 1, B n — > B and -B* — > B* strongly as n —¥ 00. Then 
for any state T £ S(J-^ N ), T Bn — > F B strongly in ©i(fj) as n — > 00. 

Proof. By Lemma [3TTl we have at least — ^ T B geometrically. However, since [r^]' 1 ' = 
(B n )[r]W(B n )* -> B[T]WB* = [T B ]W strongly in ©1(f) 1 ), the convergence of T n must be 
strong by Lemma [2T3l □ 
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3.2. Application: HVZ theorem in the general case 



In Section 12.21 we have proved the celebrated HVZ Theorem for systems with a repulsive 
interaction, W ^ 0, using the lower semi-continuity of the energy with respect to geometric 
convergence. In particular it was essential that in the absence of external field, V = 0, the 
ground state energy of the system vanishes: E°(N) = 0. 

When W has no sign a priori, the energy T i— > £ (T) is not necessarily lower semi- 
continuous, which can be seen by the fact that it may hold E°(N) < 0. Particles running 
off to infinity can carry a negative energy and in the HVZ theorem it is then necessary to 
take into account the energy of these particles. Separating the particles escaping to infinity 
from those which are bound by the external potential V is then done via localization. 

Let us recall the iV-body Hamiltonian H v (N) defined in ()2.6[) . The bottom of its spec- 
trum and of its essential spectrum are respectively denoted by E v (N) and Y, v (N). As usual 
we make the assumption that W is even, and that V and W can both be written in the form 
SiLi fi with ft G L Pi (R d ) where max(l,<i/2) < pi < oo or pi — oo but — > at infinity. 

The result in the general case is the following. 

Theorem 3.1 (HVZ in the general case). Under the previous assumptions on V and 
W, we have 



Y, V (N) =M{E v (N-k) + E°(k), k= 1,...,N}. 



(3.7) 



We now provide the proof of Theorem 13.11 This serves as an illustration of the concepts 
of geometric convergence and localization that we have introduced, but also introduces the 
reader to the techniques that we use later for nonlinear systems. 

Proof. As in the proof of Theorem 12. 1[ we only explain the lower bound ^. We take 
the same sing ular Weyl sequence {*„} such that (H V (N) - T, v (N))^ n ->■ and let T n = 
0©- • •© l^n) (^nl € J-^ N . We assume (up to extraction of a subsequence and by Lemma [2~2j) 
that r„ ^ g r = Goo © • • • © Gnn geometrically. Recall that Gnn is the weak limit of 
|^ , n)(*n|j hence Gnn = since ^„ — by assumption. 
Our goal is to prove the following fundamental estimate 

N 

S y (A) = lim E v (Y n ) > { EV ( N ~ k ) + E°(k))Tr^- k (G N ^ kN ^ k ) . (3.8) 

fc=i 

Compared to ()2.14|) . the bound now includes the energy E°(k) of particles running 
off to infinity, which can be nonzero. Recall that T is a state, that is Gkk ^ and 
X)fc=i Tr^w-* (Gjv-fc7V-fe) = 1 since Gnn — 0. Therefore the right hand side of (|3.8|) is 
a convex combination and we have 



(E V (N -k) + E (k)j r Tr S jN-k (G N - kN - k ) ^ inf {E v (N - k) + E°(k), k = 1, N} , 



N 

E 

fe=i 

which proves the lower bound in (|3.7|) . 

In order to show the inequality (|3.8[) , we pick a smooth cut-off function ^ \ ^ 1 which 
equals 1 on the ball B(0, 1) and outside the ball B(0, 2), and let xr( x ) = xi x /K) as well 
as r]n = a/1 — Xr- The rest of the proof goes as follows: 

(i) We geometrically localize in and outside the ball of radius R by means of the smooth 
partition of unity x R + Vr = 1 j 
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(ii) We use the fundamental equality (|3.6|) : 
(Hi) We pass to the limit as n — > oo; 
(iv) We take the limit R — >• oo. 

As we will explain later in the proof of Theorem 15.11 it is possible to use an n-dependent 
radius of localization R n — » oo, and to perform the steps (Hi) and (iv) simultaneously. As 
we do not need this technique here, we defer its use to Section El for pedagogical purposes. 
The so-called IMS formula reads: 

- A = xr(-^)xr + VR.(-A) m - \V X r\ 2 - |Vt7h| 2 . (3.9) 

Hence —A > xr(~^)xr + VR^^VR ~ C/R 2 . Using this for the kinetic energy as well as 
the partition of unity 1 = Xr + Vr m ^ ne interaction energy, we deduce that 



£ v (T n )^£ v ({T n ) XR )+£°((T n ) nR ) + / m (x) 2 V(x)pr n (x)dx 

W(x - y)xR{x) 2 m{y) 2 [r„] (2) (x, y; x, y) dx dy - CN/R 2 



Let us start by estimating the error terms. Since {^ n } is a Weyl sequence it is bounded 
in H 1 ((R d ) N ), thus ^/pf~ is bounded in H 1 ^), by dflTI) . Since V = J^f=i Vj with Vj € 
jj>i where max(l, d/2) < pj < oo or pj = oo but Vj — > at infinity, we have by Holder's 
and Sobolev's inequalities 

/ m(x) 2 V(x)pr n (x)dx < C V WVjrjjiW PJ( d 

which tends to zero as R — > oo. For the interaction term, we may write for instance 
W(x - y)xR{x) 2 r]B.(y) 2 [^n] (2 Hx,y;x,y) dxdy 

W(x - y)xR(x) 2 mR (y) 2 [T n ] {2 \x,y-x,y)dxdy (3.10) 

W(x - y)xR(x) 2 r lR (y) 2 X3R(y) 2 {r n ] {2) (x, y; x, y) dxdy. (3.11) 



In the first term of the right hand side, the integrand is zero except when \x — y\ ^ R, hence 
it may be estimated similarly as before by 

K 

Krai<c^ii^i w ^n^ (R<J) 

which also tends to zero when R — > oo. Summarizing we have shown that 

£ v (r n )^£ v ((T n ) XR )+£°((T n ) m ) 

W(x - y)xR(x) 2 r 1 R(y) 2 X3R(y) 2 [r n ] i2) (x, y; x, y) dxdy + e R (3.12) 



where e R is independent of n and tends to zero as R — > oo. The total energy of the system 
can be estimated from below by the sum of the energies of the localized states in and outside 
the ball of radius R, plus error terms. 
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We now deal with the main two terms and write that 

N N 

£ v ((T n ) XR )+£°((T n ) m ) =£>a* (H v (k)G« Rtk ) (ff°(fc)G™ Bifc ) 

fe=0 fc=0 

> Y J E V {k)Tr nk + J>°(fc)TV (G« R)/b ) , (3.13) 

fc=0 fc=0 

where (T n ) XR — G™ fli0 ® ' ' ' ® ^xr,n an< ^ w ^ n a similar definition for G™ Rk . At this point 
we use the fundamental relation (13. 6[) (valid since T„ is an iV-body state for all n), which 
tells us that 

Tifik (G XRk ) = Tr^w-fc (G™ RjN _ k ) 
for all fc = 0, N. Inserting in (13. 13)) and changing k into N — k in the first sum we get 

N 

S V ((r n ) XR )+£°((T n ) riR ) > Yl (E v (N-k) + E°(k))T^ N - k (G n XR ^ N _ k ) . 

k=0 

By Lemma \3. 21 (or more precisely Example 13.61) . we have (r„) XR — > T XR strongly, therefore 
lim Tr^jv-k (G x N _ k ) = Tr^ N - k {G XR , N - k ) 

where T XR = G XR , ®- ■ -®G XR , N . Recall G NN = hence G XR , N = (xr)® N G nn ( X r)® n = 
also. As a consequence, 

Af 

J^oE ( EV ( N - k ) + E°(k))Tr^- k (Gl R , N _ k ) 
k=a 

N 

= (e v (N -k) + E°(k))Tr^- k (G XR , N . k ) . 
fe=i 

Using that the term in (I3.1ip converges as n — > oo since XR( x ) 2r lB.(y) 2 X3R,(y) 2 nas a compact 
support, we arrive at the estimate 

N 

H V (N) = lim £ v (T n ) ^ V (E v (N-k) + E°(k))TY^ N - k {G XR , N . k ) 

+ 2 [ [ W(x-y)xR(x) 2 r) R (y) 2 x 3R (y) 2 [r}W(x,y;x,y)dxdy + e R . (3.14) 

Passing finally to the limit R oo (using that T XR — > T strongly by Lemma 13.31 hence 
G XR} k — ► Gfcfc as i? — > oo) gives the desired estimate (|3.8p and ends the proof. □ 

4. Finite rank approximation of many-body systems 

In the previous two sections we have introduced geometric tools for many-body systems and 
we have illustrated their use on linear systems (the HVZ theorem). In practice, physicists 
and chemists resort to approximate models which are simpler to handle and to simulate 
numerically. These approximations are usually classified in two different categories: 

• those in which the set of states is reduced, 

• those in which the energy is modified by adding nonlinear empirical terms. 
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These two methods can of course be combined: in the so-called Kohn-Sham method of atoms 
and molecules 28 , all states are assumed to be of Hartree-Fock type but the energy is further 
modified to take into account exchange-correlation effects. Both techniques usually lead to 
nonlinear models, either because the class of states is replaced by a manifold or because the 
energy is itself nonlinear. 

The purpose of this section is to study methods of the first kind in which the many-body 
energy is kept linear, but the set of states is reduced. Methods from the second category 
will be considered in Section [5j We study here the so-called finite rank approximation which 
consists in assuming that the TV-body wavefunction can be expanded as tensor products 
of finitely many unknown one-body functions {ipi, <p r }. For fermions, this leads to the 
celebrated Hartree-Fock method 35J when r — N, and to the widely used multiconfiguration 
methods [19, 30] when r > N. For bosons, the Hartree method is obtained when r = 1. 

We investigate properties of geometric limits of finite-rank states, and deduce nonlinear 
versions of the HVZ Theorem. As we will see, the situation is however still rather unclear for 
bosons and our results are only satisfactory for fermions in the Hartree-Fock approximation 
or for multiconfiguration methods with repulsive interactions. We hope to come back to the 
other interesting cases in the future. 

4.1. States living on a subspace of Sj, finite rank states 
4.1.1. Definitions 

Definition 4.1 (States living on a subspace). Let Sj' C fj be a closed subspace of the 
one-body space f) 1 and P be the orthogonal projection onto Sj'. A state Y g S(J-^ N ) is said 
to live on S)' when Yp = Y. 

The smallest subspace ij' such that Y lives on US,' can be called the support of Y. The 
following is a reformulation of a result of Ldwdin |43j stating that the support can be found 
by means of the one-body density matrix [T]W only. 

Lemma 4.1 (Lowdin's criterion). Let Y be a state on J 7 ^ 1 ^ and P : f) 1 — > f) 1 be an 

orthogonal projector. The following assertions are equivalent: 

(1) Y lives on Pfy , that isYp = Y; 

(2) p[r]«p= [r]« ; 

(3) PrP = r where F = 1 ® P ® (P <Z> P) ® ■ ■ ■ ® P® N . 

Proof. It is clear from the definition of geometric localization that (1) implies (2). If we 
denote by G k i the matrix elements of Y, (3) means that P® k G k tP® = Gki for every 
^ k,£ ^ N. Using (Ti~T9l) , this is easily seen to imply that P»p\T]M P®i = [Y]^ for all 
^ p, q ^ N, hence Yp = Y and (1) holds true. 

It therefore only remains to show that (2) implies (3). We denote as usual by Gki the 
matrix elements of Y and note that, by (|1.19|1 . 

N 

[Y]M = J2[G kk ]W. 

k=l 

Our assumption that P[Y]^P = [T]W implies that 

P[G kk ]Wp= [G kk ]^ for allfc = l,..,/V. (4.1) 
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Indeed, we have P^jWp 1 = = J2 k =i P X [Gkk] W '-P X where P x = 1-P. Since [G kk ]W ;> 
for all k = l,...,N, this implies that P A -[G k k] {1) P X = 0. Now dHJ) follows for instance 
from the fact that 

P X [G kk ]^) (p x [G fefe ]«)* < P^G^ 1 ^ = (4.2) 

which shows that P x [G kk ]W = [G k k] W P X = 0. 

We now prove that flO) implies that P® k G kk P® k = G kk . We have for any P 2 , P fe G 

TV (P x ® P 2 ® ■ ■ ■ P k G kk P x <E> P 2 ® ■ ■ ■ P fe ) 

< Tr^ (P x <g> 1 ® •• • lG kk P x ® 1 ® •• -1) = -Tr^ (p^IGfefe]^)) = 0, 

by (|1.15|) . The argument is the same if P 1 - is not in the first place of the tensor product but 
appears at another position. Arguing as before, this implies P® k G kk P® k = G kk . For the 
off-diagonal terms, we have G k gGi k ^ G kk since ^ T ^ 1. This can be used to show that 
P x ® P2 ® • • • ft G*/ = 0, hence P® k G u P® 1 = G u . □ 

We now use the previous concept to define finite-rank states. 

Definition 4.2 (Finite rank states). A state T G S(J-^ N ) is said to have a finite rank 
when it lives on a subspace of finite dimension, that is when there exists a projector P of 
finite rank such that Tp = T. The rank ofT is then defined as 

rank(r) = min{rank(P) : T P = T, P 2 = P = P*} = rank (V] (1) ) . 

The last equality follows from Lemma T4.1 I 



Example 4.1 (Coherent, Hartree and Hartree-Fock states). For bosons, both the 
Hartree state \(p® N ) and the coherent state W(f)\SY) have rank r — 1. For fermions, a pure 
Hartree-Fock state (p± A • • • A tpjf has rank r = N. o 

The following says that finite-rank states are dense in S(J-^ N ). 

Lemma 4.2 (Approximation by finite-rank states). Any state T G S(J-^ N ) is a strong 
limit of finite rank states. 

Proof. Let {<fj} be an orthonormal basis of fj and P„ := 53j=i If^ifjl- Then P n — » 1 
strongly in $). Therefore by Lemma I3~3l it holds Tp n —> T strongly. But Tp n has finite rank 
since (Tp n ) p — ^(p n y — Tp n by Remark l3.2l □ 

We now show that any state of finite rank is a finite linear combination of monomials in 
the creation and annihilation operators. 

Lemma 4.3 (Expansion of finite rank states). Assume that Tp = T for some or- 
thogonal projector P = Y^j=i IVjX^il of finite rank r, and let (G k e)i<^ ki ^N be the matrix 
elements ofT. Then each G k g can be expanded as follows: 

G kl = c u a\L Pll )---a){L Plk )\Vl){Vl\a(ip H )---a{ip n ) (4.3) 

I={*i^-<»fc}C{l,...,r} 
J={ji<-'-<Ji}C{l,...,r} 

for some Cij £ C. 
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Proof. This follows from the fact that P m G M P® e = G M , see (3) in LemmaO □ 

Consider a finite rank state, that is such that [r]*- 1 ' has finite rank r (Lemma l4.ip . Then 
we can write [r]^ 1 ) = Y^j=i n j\ i Pj)( l Pj\ f° r an orthonormal system {<fj}j—i of eigenvectors of 
[p]( 1 ) > xhe rij are usually called the occupation numbers and the ipj the natural orbitals of 
r. Lemma l4.3l then shows that any finite rank state can be expanded by means of its natural 
orbitals. This is the original version of Lowdin's Expansion Theorem [43] (see also Lemma 
1.1 (ii) in [19] and Lemma 1 in [30]). 

The simplest example is that of a state of the form T = © • • • © |^'}(^ r |, that is a pure 
iV-body state. Then if rank([r]( 1 ^ ) ) ^ r and {ipi, <p r } is an associated orthonormal system 
of natural orbitals, it holds 

* = ^2 cn,...,i N <Ph ° ■ ■ ■ ° <Pi N 

where o = A (fcrmions) or V (bosons). 

4.1.2. Geometric properties of finite rank states 

We now turn to properties of finite-rank states with regard to geometric localization and 
convergence. The following is a simple consequence of the characterization of the rank in 
terms of the one-body density matrix (Lemma 14. ip . 

Lemma 4.4 (Localization and geometric limit of finite rank states). 

(1) If a state T G S(J-^ N ) has rank ^ r, then for every localization operator B, ^ 
BB* 1, the corresponding localized state Tb has rank ^ r. 

(2) If {r„} is a sequence of states on of rank ^ r and T n —*- g T geometrically, then 
r has rank ^ r. 

Proof. The result follows from the fact that when rank([r] (1) ) r, then rank (B[T] (1) B*) ^ 
r for every localization operator B. Similarly, when [Tn]^ — s [r]^ 1 * 1 weakly-* in $), then 
rank([r]«) s$ liminL^ rank^^ 1 )). □ 

For iV-body systems we often have to study sequences of states of the form T„ = ffi 
• • • © l^nX^nl- When r„ —*- g r geometrically and when each T n has rank ^ r, then we have 
by Lemma l4~4l T = Goo © • • • © Gnn where each Gkk has rank < r. A similar property holds 
for a localized state T b ■ This information is unfortunately not enough to be really useful in 
applications. It is fortunate that this can be precised in the fermionic case, as expressed in 
the following important result. 

Lemma 4.5 (Localization of a fermionic TV-body finite-rank state). Let G <E S(^) 
be a fermionic state of the N -body space Sj^ 1 , of rank ^ r, and r = 0ffi---©0©G be the 
corresponding state in J 7 ^ ■ Let B be a localization operator, ^ BB* ^ 1, and denote by 
Tb = Gq © • • • © G^n the corresponding localized state in "F^ ■ Then each G^ k belongs to 
the convex hull of k-body states of rank at most r — N + k: we have 

G kk = ^2 a> i S 3 

3 

with Sf € S(Sj-%), ctf^O and 

rank(S*) < r - N + k. 
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This result does not hold in general for bosons. In Example 13.41 we have seen that the 
localization r = Gq © • • • © G% N of a Hartree state r = © • • • © \ip® N ){<p® N \ with rank 
r = 1 satisfies rank(Gf fc ) = 1 for all k = 1, N. We now provide the proof of Lemma \A. 5 1 



Proof. Since G has rank at most r, there exists a projector P = Y]j—\ [tp^ifj] oi rank r 
such that Tp = T. By linearity we can assume that G is a pure state, that is G = | V E , )(5'| 
where 

*= ^2 Cii-iN<Pii<8>---®<Piir- (4-4) 

We follow here the notation of [30] : <ki...i N reflects the symmetry of the wavefunction, that 

is c i CT( i)...i CT (jv) = e ( cr ) c ii...iN an d Cjj ^ = as soon as two indices are equal. 

We have the freedom to choose any orthonormal basis of the finite-dimensional space 
V = span(</3j) = Range(P). Indeed, if we replace the functions <pj by <p'j — J2i-iUij<pj 
for an r x r unitary matrix U = (Uij), then l|4.4p is still valid, with adequately modified 
configuration coefficients c' i± ijv (see Formula (12) in [30]). Taking advantage of this gauge 
freedom, we can diagonalize any matrix of the form (tpi,Mipj) for a well-chosen one-body 
self-adjoint operator M. Here we choose M = B*B, that is we work with an orthonormal 
system for which the r x r hermitian matrix ((Bipi, Bipj))±^i j<^ r is diagonal. In particular 
we have (Btpj,Btpj) = Q if i ^ j, which dramatically simplifies the expression of Tg. Using 
Formula (|3.4jl . we find that 



G% k =( I Pi y fi-iiB^ +l f)--(i -i^.vini*... ,...!*><*.. ,..,1 

where 



Ve k+1 ...e N = c ii-i k e k+1 ...e N Bip n ® ■ ■ ■ ® Btp ik 

ill-- ,ik&{l,---,r} 

y Ci 1 ...i h t h+1 ...e. N B(p h ®---®B<p ih . (4.6) 

il,— ,ik€{l,...,r}\{l k+1 ,...,e N } 

In (|4.6[) . we have used that, for fermions, i N = when two indices coincide. Clearly, 
^t k+1 ...i N has rank ^ r — N + k and the result follows. □ 

Example 4.2 (Localization of pure Hartree- Fock states). Let ^> := <p\ A • • • A <^jv be 

a pure Hartree-Fock state and T := © • • • © | v I')(vE , | be the corresponding state in J 7 ^ ■ 
The localization of T is T B = Gfi © • • • © G% N , with 

/ i< **,}< 1 1 v: 

JJ (l-|^W) A-A%)(% A---AB<p ih \, (4.7) 

e{i,...,N}\i J 

assuming that the orbitals have been chosen such as to ensure {Bipi, Bipj) = when i ^ j.o 

From Lemma 14.51 we can deduce the general form of the geometric limit of fermionic 
A^-body finite-rank states. 

Lemma 4.6 (Geometric limit of fermionic A^-body finite-rank states). Let T„ = 

© • • • © G n G J-^ N be a sequence of fermionic N -body states, with rank(r„) r for all n. 



Geometric methods for nonlinear many-body quantum systems 



35 



If r„ — ^ g r = Goo © • • • © Gnn geometrically, then each Gkk belongs to the convex hull of 
k-body states of rank at most r — N + k. 



Proof. Let {(fi} be any fixed orthonormal basis of $j and Pj := Ylj=i IVjXwl be the 
projector onto the space spanned by the first J elements of this basis. Since Pj is compact 
for every fixed J, we have by Corollary 13.11 (T n )pj — > (T)pj strongly. By Lemma |4"31 we 
know that each (L„)p 7 can be written in the form (L„)p ; = ®^L G fc j^, each G^™ being a 
convex combinations of states of rank at most r — N + k. By strong convergence, we infer 
that (r)p 7 has the same property. Now since (r)p J — > L strongly as J — > oo by Lemma 1331 
we conclude that L also satisfies the same property. □ 



Example 4.3 (Geometric limit of pure Hartree-Fock states). Let <F n := <y5™A- • -A<p N 
be a pure Hartree-Fock state and L„ :— 0®- • - ® | x I'n)(vE , n| be the corresponding state in . 
We assume that — ifj weakly in f), for j = l,...,N. Up to applying an n- independent 
unitary transform U to the <£>™'s we may also suppose that (<pi, tpj) = when i =/= j. We then 
have that L„ —*- g Goo (S ■ ■ ■ (3 Gnn geometrically, with 

Gkk = Yl n (i-IIVal 2 ) ] l^ii A---Ai/3i fc )(v? il A---Atp ik \. 

Z={i 1 <-<i fe }C{l,...,iV} \a£{l N}\I J 

(4.8) 

We see that 

either there is strong convergence, ip™ — > ifj for all j — 1, ...,JV, hence Gkk = for all 
k = 0, ...,N- 1; 

or all the particle are lost, <pj = for all j = 1, N, thus Gkk = for all k = 1, N, that 
is L = |0)<0|; 

or not all the particle are lost, that is < \\ipj\\ < 1 for at least one ipj, and there exists 
1 k ^ N - 1 such that G fefe ^ 0. 

Indeed if we assume (up to reordering) that ipi, ...,ipN 1 ^ but <PNi+i = ■ ■ ■ = lpn = 0, 
we see that Tr{GN 1 N 1 ) ^ Y\.j=i ll^j'll 2 > 0- The fact that we cannot have Gkk = for all 
k = 1, N — 1 while both Goo and Gnn are ^ will be very useful later in the proof of 
Theorem W5\ o 



4.2. HVZ-type results for finite-rank many-body systems 

4.2.1. A general result 

Let us come back to the TV-body Hamiltonian 

N / A \ 
H V (N) = J2(—^ + V( Xj ))+ Y, W(x k -xt) 



which we have already introduced in (|2.6|) . As usual we make the assumption that W is 
even, and that V and W can both be written in the form J^Lj fi with /j S L Pi (IR rf ) where 
max(l, d/2) < p i < oo or = oo but /j — > at infinity. 

For bosons or fermions we may introduce the approximated ground state energy obtained 
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by restricting to finite-rank states: 

EV(N):= inf (*, H V {N)V). (4.9) 

*6 ff a/»(( R ) ) 

rank(*)^r 
11*11=1 

We clearly have E?{N) ^ E v (N) for all r, and linWoc E^[N) = E V (N). 

Let us emphasize that, although the energy functional is the same as for the full linear 
model, we now have the additional constraint that rank(^) ^ r which is itself highly nonlin- 
ear. Thus the so-obtained Euler-Lagrange equations are themselves nonlinear. If r = 1, one 
gets for bosons the Hartree nonlinear equation. For fermions, one obtains the Hartree-Fock 
equations [35l [42] for r = N and the multiconfiguration equations [HI [30] for r > N. 

We are interested here in existence results for ground states by means of geometric meth- 
ods. The following theorem is a generalization to the nonlinear case of the HVZ Theorcm l3.ll 

Theorem 4.1 (Finite rank HVZ-type result, general case). If the following inequal- 
ities hold true 

EV(N)<E?(N-k) + E°(k), Vk = l,...,N, (4.10) 

then all the minimizing sequences {^ n } for the variational problem E^f (N) are precompact, 
hence converge, up to a subsequence, to a ground state of rank r. 
If all the particles are fermions, (|4.10j) can be replaced by 

E^(N)<E^_ k (N-k)+E r _ N+k (k), Vfc = l,...,7V. (4.11) 

Proof. Up to a subsequence we may assume that * n * weakly in H^,((M. d ) N ), and 
that the corresponding state T n := © • • • © |4 r 71 )(^ r n | € S(J-^ N ) converges geometrically 
to T = Goo © • • • © Gnn- If ||^[ = Tr(Gjvw) = 1 then we have strong convergence r„ — s- V 
in &i(T^ N ), hence $ n — > W in I? . Under our assumptions on W, this can then be used to 
prove that the two-body term converges strongly: 

lim / dxi---f dx N W(xi - Xj)\^> n (x!, ...,x N )\ 2 

= V] / dxi ■ ■ ■ I dx N W(xi - Xj)\^>(xx, ...,x N )\ 2 . 

Since the interaction term is the only one which can fail from being weakly lower semi- 
continuous, we deduce that 



E?(N) = lim £ v (* n ) > £ v ^>) ^ EY(N), 



hence that ^ is a ground state for E V (N). Finally, strong convergence in if*, ((R d ) ) is 
obtained by noting that lim JWOO £ v ($ n ) = £ v (^f), hence that the kinetic energy must also 
converge. 

Summarizing the previous paragraph, we only have to prove that Gkk — for all k = 
0, N — 1. We follow the proof of Theorem 13.11 we localize the system in and outside a 
ball of radius R, by means of a smooth partition of unity, \r + = 1- In the lower bound 
corresponding to (|3.13[) , we may use that each G™^ k has rank ^ r by Lemma 14.41 (or rank 
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^ r — N + k for fermions, by Lemma l4~5]) . To be more precise, each G" R k can be diagonalized 
as follows 

/~m \ R,k,n i T R,k,n\ i T R.k.n i 

i 

where g^' k ' n ^ and (p n )® fc i[/^' fc '" — f or an orthogonal projector P n of rank ^ r 

(or r + N — k for fermions) . Saying differently each G™ fl k is a convex combination of pure 
states of rank ^ r. Hence we have an estimate of the form 

TV (H V (k)G^ k ) > E^{k) TV (<w) - 

with Ey (k) replaced by E^_ N+k (k) for fermions. A similar argument applies to the terms 
involving G" fl fe . Taking the limit n — > oo first and then removing the radius R of the 
localization, following the proof of Theorem 13. 1[ we arrive at the following estimate, similar 
to (EES)) : 

iv 

<(A0 > £ (i? r y (fc) + E° r (N - k)) Tr^(G fefe ) 

k=0 

(with an obvious modification for fermions) . The term on the right is a convex combination 
of E¥(N) (for k = N) and E?Qe) +E°(N - k) for k = 0, N - 1. When (|4T0j) holds, this 
is only possible if Gkk = for all k — 0, 2V — 1. □ 

Unfortunately Theorem 14.11 only provides a sufficient condition for the compactness of 
minimizing sequences. In general we do not expect that (|4.10l) (or (|4.11[) for fermions) is 
also a necessary condition. The reason is that when two systems are placed far away in 
space, the rank of the whole system becomes the sum of the ranks of the two subsystems. 
This sum being 2r for (|4.10[) and 2r - N for (14. lip , the inequalities (14.10)) and (|4.11[) are 
not expected to be correct in general when the strict inequality < is replaced by a large 
inequality ^. It is usually when large inequalities hold true that one can get necessary and 
sufficient conditions. 

In the next section we will give two examples for fermions, due to Friesecke [19] , for 
which one can reduce (|4.11l) to inequalities of the form 

E?(N) < E^_ r ,(N - k) + E°,(k), (4.12) 

hence providing a necessary and sufficient condition of compactness of minimizing sequences. 
The case of geometric methods for finite-rank bosonic systems is still largely unexplored. 

4.2.2. Two corollaries for fermions 

We give two corollaries of Theorem 14. II in the fermionic case. These two results are contained 
in a paper [19] of Friesecke (see in particular Corollary 6.1 of [H]), with a proof that is not 
very much different from our approach. Our formalism automatically takes care of the com- 
plicated geometrical methods for finite-rank states which was detailed in [19] (in particular, 
the reader should compare Friesecke's Lemma 4.1 in !19 with our Lemma [43]). 

The first result deals with the Hartree-Fock case, corresponding to having rank r = N. 

Corollary 4.1 (Hartree-Fock HVZ-type). Assume that all the particles are fermions, 
and that V and W satisfy the same assumptions as before. Then the following assertions are 
equivalent: 
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(1) E%(N) < Ejf_ k (N - k) + E°(k) for all k=l, —,N; 

(2) all the minimizing sequences {^ n } for the Hartree-Fock ground state energy E^(N) 
are precompact in H^((M. d ) N ), hence converge, up to a subsequence, to a minimizer for 
EV(N). 

Proof. The implication (1)=> (2) follows from Theorem 14.11 in the fermionic case, with 
r = N. To prove the converse inequality we first notice that it always holds E^(N) ^ 
EYf_ k (N — k) + E k (k) for all k — 1, ...,N. This is easily seen by taking a trial function of 
the form 

*„ = I- 1 Af 2 (- -nv) (4.13) 

where v € K d \ {0}, = A ■■ ■ A (fiN-k and & 2 — tpN-k+i A ■ ■ ■ A <pn are trial functions 
for, respectively, the problems E^_ k (N — k) and E k (k). For simplicity one can take all the 
ipj's of compact support. If there is equality Ej^(N) = E^-_ k (N — k) + E k (k) for some 
k G {1, N}, then a minimizing sequence for Ej^(N) of the same form as (|4.13[) can be 
constructed and it is clearly not compact. This shows the converse implication (2)=>(1). □ 

There are now many different proofs for the existence of ground states in Hartree-Fock 
theory. For atoms and molecules, the first is due to Lieb and Simon [35] ■ An approach 
based on a second-order Palais-Smale information was proposed later by Lions [32]. These 
two methods rely on a formulation of the problem in terms of the N orbitals tpx, ...,ipj\[ 
of the Hartree-Fock state as well as on the assumption that W ^ 0. A different approach 
due to Lieb [32J (see also [2J UJ |3]) uses generalized Hartree-Fock states and the fact that, 
when W 0, a generalized ground state is necessarily a pure state. In this formulation the 
minimization problem is expressed using as main variable the one-body density matrix [r]^ 1 -* 
which completely characterizes the Hartree-Fock state. When W is not positive, it cannot be 
guaranteed that a generalized ground state is necessarily a pure state, and Lieb's variational 
principle of [32] cannot be employed. 

Our approach here (due first to Friesecke [19]) is completely different and it is based on 
geometric properties of TV-body Hartree-Fock states. It leads to quantized inequalities of the 
form of that of Corollary 14. 11 without any assumption on the sign of W. 

Of course, the next step when studying a specific model is to prove that the binding 
inequality holds true. As explained by Friesecke in [19] , this can be done by induction: using 
that there exist ground states for the problems with k particles (1 ^ k < N), one tries to 
prove by a convenient trial state that E^ (N) < E]^_ k (N — k) + E k (k), showing the existence 
of a ground state for E^(N). For atoms and molecules, this argument can be carried over 
as soon as N — 1 < Z, where Z is the total charge of the nuclei. 

Our second application of Theorem 14.11 in the fermionic case is the multiconfiguration 
case N ^ r for repulsive interactions. 

Corollary 4.2 (Multiconfigurational HVZ-type in the repulsive case). We assume 
that all the particles are fermions, that V and W satisfy the same assumptions as before 
and, additionally, that W 0. For every r is N , the following two assertions are equivalent: 



(1) EY{N)<EY^{N-l); 
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(2) all the minimizing sequences {^ n } for E^ (N) are precompact in iJ^((R ) ), hence 
converge, up to a subsequence, to a minimizer for E*. (TV) . 

The reason why we restrict to W ^ is because it then holds 

E° r -N + k(k) = E° k (k) = 0. 

Hence if we insert this in (|4. 1 1[) we are left with an inequality of the form of (|4.12[) . It is still 
an open question to understand the geometric behavior of multiconfiguration methods for 
non- repulsive interaction potentials (see, in particular, the comments on page 56 of |19]). 

Proof. The proof follows that of Corollary 14.11 using that E°_ N+k (k) = since W ^ 0, 
and that mi k=h ..., N {E^ k (N - k)} = E^(N ~ I). □ 

Again for atoms and molecules, one can prove by induction the existence of a ground 
state as soon as N < Z + 1, see [19]. 

4.2.3. Translation-invariant Hartree-Fock theory 

In this subsection we study a translation-invariant Hartree-Fock model, that is we assume 
that V = 0. It is known that (by translation- invariance) the A-body Hamiltonian H°(N) 
never has any ground state, but it can happen that there is one when restricting to Hartree- 
Fock states. Of course translation-invariance is not really broken: minimizers are not unique 
as they can be translated anywhere in space and it is the whole set of minimizers which is 
invariant under translations. 

Because of the action of the group of translations it can only be hoped to prove com- 
pactness of all minimizing sequences up to translations. 

Theorem 4.2 (Translation-invariant Hartree-Fock). We assume that W satisfies the 
same assumptions as before (but W need not be non-negative). Then for all N ^ 2, the 
following assertions are equivalent: 

(1) E° N (N) < E%_ k (N - k) + E° k (k) for all k = l, N — 1; 

(2) all the minimizing sequences {^ n } for E^(N) are precompact in H^((M. d ) N ) up to trans- 
lations. Hence there exists {v n } C M. d such that — "„) converges, up to a subsequence, 
to a Hartree-Fock minimizer for Epj (iV) . 

The notation *&„(• — v n ) is interpreted in the sense of (xi, xn) ^ ^(^l — v n , xn — 
v n ). A result of the same kind was shown for the first time by Lenzmann and the author 
in [2 9) , for a model of neutron stars with a pseudo-relativistic kinetic energy and the grav- 
itational Newton interaction. The pseudo-relativistic kinetic energy yields new difficulties 
concerning boundedness from below of the energy and localization errors (see Lemma A.l 
in [29]). For nonrelativistic systems one easily arrives at the following result: 

Corollary 4.3 (Nonrelativistic Newtonian Hartree-Fock systems). Assume that all 
the particles are fermions, thatd = 3 and W(x— y) — — g/\x— y\ with g > 0. Then E^(N) has 
a Hartree-Fock ground state for all N ^ 2 (hence infinitely many by translation-invariance) . 

Proof. The binding inequality Epj(N) < Ef f _ k (N — k) +E k (k) can be proved by induction 
using Newton's theorem, as explained in 29 . □ 
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We are now ready to prove Theorem 14.21 

Proof. It was already shown in the proof of Corollary 14.11 that E° N (N) ^ E^_ k (N — k) + 
E®(k) for all k = 1, ...,N — 1. Furthermore, if there is equality for some k, then one can 
construct a minimizing sequence which is not compact, even up to translations. Therefore 
we only have to prove that (1)=>(2). 

To this end, we consider one minimizing sequence ^ n = (p n A • • • A ip 1 ^ for (N) and 
we define the associated state in , T n = © • • • © |4 r n )(^„|. Since £ °(r„) is bounded, 
by (|2.7p we have a uniform bound on the kinetic energy: 

Tr4(-A)[r„p)) <C. 

This itself implies a uniform bound on the H 1 (R d ) norm of y/pr n , by the Hoffmann-Ostenhof 
inequality (|2.1ip . 

Our goal is to prove convergence of 4 r n (- — v n ) for an appropriate translation v n . The first 
step is to determine this translation v n by detecting a piece of mass which retains its shape 
for n large and, possibly, escapes to infinity. We therefore consider all the possible geometric 
limits, up to translations, of subsequences of {T n } and we define the largest possible average 
particle number that these limits can have: 

m({r„}) := sup jTr^(Afr) : 3{v k } C R d , r^r„,r_^ T J. (4.14) 

Here t$ is the translation unitary operator defined by (t^)(xi, xn) = ^(xi—v, xn~v) 
when \1/ £ Sj N and extended by linearity on the whole Fock space. By the strong convergence 
pr„ — > Pt in ^ioc(^ d ) when r„ -± g T (with bounded kinetic energy), we also have that 

m({r„}) = sup | p : 3{t> fe }cK d , p r „ k (• - « fe ) 1/2 P 1/2 weakly in H 1 (M d )|- (4-15) 

The definition of m({r n }) is inspired of a result of Lieb [33] as well as of the concentration- 
compactness method of Lions [38l [39] . The purpose of m({T n }) is to detect the piece contain- 
ing the largest average number of particles, which possibly escape to infinity (when \vk\ — > 
oo). Following Lions' terminology, a sequence {F n } is said to vanish when m({T n }) = 0, 
which is equivalent to the property that 

g 

or that 

V{4}ci d , pr n {- - v n ) ^ &.e. 

As we now explain, saying that m({r„}) = is actually quite a strong statement. 

Lemma 4.7 (Vanishing). Let {T n } be any sequence of states on J-^ N , with a uniformly 
bounded kinetic energy. The following assertions are equivalent: 

(i) m({r„}) = 0; 

(ii) for all R > 0, one has lim sup / pr n = 0; 

n ^°°xm d JB(x,R) 

(Hi) pr n strongly in L p (R d ) for all 1 < p < p* , where p* = d/[d— 2) if d 3, p* = oo 
if d= 1,2. 
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Proof. The fact that (i) => (ii) follows from the strong local convergence of pr„ ■ The 
implication (ii) =>• (Hi) was proved first by Lions in }39j (Lemma 1.1). Finally, it is clear 
that if pr n — > strongly in one L p (M. d ), then pr n (' — x n ) — ► strongly in L p (M. d ) for every 
sequence {x n } C R d , hence (i) follows. □ 

We will now show using Lemma 14.71 that our Hartree-Fock minimizing sequence cannot 
vanish. We have, using Wick's Theorem for generalized Hartree-Fock states [351 II] > 

*n, I E W(x i -X j )\^ n 

<(*», E l^l^i-^))^ 



= W / l^-y)l(pr„(ik„(3/)-|[T n ] (I) ^w)| 2 ) ^ 
< \ [ Pv n (pr n * \W\) . 



When m({r„}) = 0, we have that p r „ ->• in L p (R d ) for all 1 < p < p* by Lemma l4~7l 
Under our assumptions on W, this implies that the interaction term converges to 0. The 
kinetic energy being non-negative, this shows that in the case of vanishing 

E° N (N) = lim £°(T n ) > 0, 

n— >oo 

which contradicts the assumption that (1) holds true (it is clear that (1) implies that 
E%(k) < kEl(l) = hence, since the inequality is strict in (1), that E%(N) < 0). 

We have shown that m({r„}) > 0. This proves that there exists a sequence {vk} C R d 
and a subsequence r„ fc such that T' k := T.$ k r nk T_z k — ^ T with T ^ |f2)(f2|. Since the 
problem E^(N) is invariant under translations, the new sequence T' k is also a minimizing 
sequence for E^(N). To simplify our exposition, we do not change our original notation and 
we assume that T n ^ g T with T — Goo©' ■ '©Gnat- The assumption that T =/= means 
that ^ Goo < 1. As usual strong convergence of {^ n } in L 2 implies strong convergence in 
H 1 and it suffices to prove that Gkk = for all k = 0, N — 1. 

We can now follow the proof of Theorem 14.11 which uses a localization in a ball of radius 
R as well as strong convergence in this ball, before passing to the limit R — > oo. This yields 
an inequality of the form 

N 

E° N (N) > J2 ( E °N-k(N - k) + E° k (k)) Tr fl «- fc (G fcfc ). 

k=0 

Note that in comparison with Theorem 14.11 the terms corresponding to k = and k = N 
are equal. When the binding inequality holds, this is only possible when Gkk = for all 
k = 1, ...,N — 1 . Hence we have 



r = G 00 ffiO©---©o©G 



NN- 



We also know that Goo 7^ L hence Gnn 0. We have already explained in Example! 
that the only geometric limit of a sequence of pure Hartree-Fock states of this form must 
have Gqo = 0. This ends the proof of Theorem 14.21 □ 
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5. Many-body systems with effective nonlinear interactions 

In this section we consider a system of N quantum particles whose many-body energy is not 
linear with respect to the state | X I')( V I / | of the system, but also contains a nonlinear term F: 

= (*,fr(JV)*) + F(|*)<*|). 

The purpose of the last term is often to effectively describe complicated interactions between 
our TV particles, through a second quantum system which has been eliminated from the 
model. Even when the model is translation-invariant, the N particles can form bound systems 
thanks to the nonlinear term F. 

Situations of this kind are ubiquitous in quantum physics. In Section l5.2[ we study the 
example of the ./V-polaron, which is a system of N electrons in a polar crystal. In the so-called 
Pekar-Tomasevich model, the crystal is eliminated and replaced by an effective nonlinear 
Coulomb-like force between the electrons. 

In nuclear physics, strong forces between nucleons are also often described by effective 
nonlinear terms. The most celebrated ones are the Skyrme [S3] and the Gogny fTU] forces. 
Although these methods have been mainly used in the context of mean-field theory, their 
extension to correlated models was recently considered in [35] . 

In this section we illustrate our geometric techniques by studying the simple case of 
a concave nonlinear term F depending only on the density py of the system. We state a 
general theorem in Section 15.11 and apply it to the multi-polaron in Section 15.21 



5.1. A general result 

Let us consider a system of N spinless particles (bosons or fermions) in E rf , interacting via a 
potential W and a nonlinear effective term F. For simplicity we assume that F only depends 
on the density of charge p^ of the many-body state 




-A, 



^2 W{x k - xt) 



We also introduce the corresponding ground state energy, for bosons or fermions, 

E(N) = inf £(*). 

*Gff^ /s ((R d ) JV ) 
11*11 = 1 



(5.1) 



(5.2) 



As before we make the assumption that W can be written in the form Ylif=i W< with 
Wi £ L Pi (R d ) where max(l,d/2) < pi < oo or pi = oo but W% at infinity. As for the 
functional F, we assume that it satisfies the following assumptions: 

(Al) (Subcriticality) F is a locally uniformly continuous functional on L Pl (R d ) n L P2 (R d ), 
for some 1 < p% ^ pi < p* , where p* = d/(d— 2) when d > 2 and p* = oo when d = 1, 2, 
and such that F(0) = 0. Furthermore, there exists < e < 1 and C > such that 



Vip £ H\R d ), 



M 2 ^ 



%| J )>-5 / |V^| 2 -C 



(5.3) 



(A2) (Translation invariance) F(p(v 
v e R d ; 



■)) = F(p) for all p G W^W 1 ) n L P2 (M d ) and all 
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(A3) (Decoupling at infinity) If {pl^} and {p^} are two bounded sequences of L 1 (K d ) n 
£P 2 (R d ) such that d(supp(pi) , supp (/>*)) — > co, then it holds 



(A4) ^Concavityj F is concave on the cone {p G i Pl (M d ) n L P2 (R d ) : p ^ 0}; 

(A5) (Strict concavity at the origin) For all p € L Pl (K d ) n L P2 (M d ) with p > and p ^ 0, 
one has F(ip) > for all < t < 1. 

Example 5.1. Consider the following functional: 



It can be verified that F satisfies all the previous assumptions when a>0,l</3<l + 2/<i, 
and when the function h is of positive type (h ^ 0) and can be written in the form h — 
Si=i hi with hi £ L qi (R d ) for some max (l, (d + l)/2) < qi < co. When d = 3, this covers 
Coulomb interactions h(x) — l/|x|, as well as Dirac's term corresponding to (3 = 4/3. o 

In the proof, the concavity of the functional F is crucially used to extend the energy 
£ to mixed states in the truncated Fock space , making possible the use of geometric 
methods. Concavity might seem a very strong assumption but it is indeed very natural from 
a physical point of view. As we have explained the term F(p^) usually empirically describes 
the interaction of our N particles with a second (infinite) system (for instance phonons of 
a crystal for the multi-polaron studied in Section |5.2[) . In most physical models the real 
coupling between the two systems is linear with respect to the state of the N particles (for 
instance linear with respect to py). Eliminating the degrees of freedom of the second system 
by simple perturbation theory or minimization over product states always leads to concave 
functionals F. 

The assumption (Al) that F is subcritical will be used in the proof to discard the possi- 
bility that minimizing sequences vanish. The other assumptions on F are of a more technical 
nature, and they can certainly be relaxed a bit. It is possible to treat non translation- 
invariant functionals but in this case the main result below is not stated the same. It is 
also easy to generalize the main theorem below to the case of a functional F which is not 
a simple function of the density (for instance when F is a function of the one-body density 
matrix), with appropriate assumptions. 

It is a simple exercise to verify that, under the previous assumptions, the energy func- 
tional £ is well-defined and continuous on H*, a ((Mr) N ). Moreover, using (|5.3|) in (Al) and 
the Hoffmann-Ostenhof inequality (|2.11l) . we have 



F(Pn + Pi) ~ F(Pn) - F(pl) -> as n -> co; 





(5.4) 



In the second line we have used the assumptions on W, similarly as in (|2.7|) . This shows 
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that £ is bounded from below, hence that E(N) is finite. In the following we denote by 



N 



-A 



H(N) := ]T + J2 W ( x * ~ x ^ 



the translation-invariant many-body Hamiltonian. The main theorem is the following: 

Theorem 5.1 (Nonlinear HVZ for many-body systems). Under the previous assump- 
tions, the following assertions are equivalent: 



(1) One has 



and 



E(N) < E(N -k) + E(k) for all k = 1, N -I, 



E(N) < m£a(H(N)) : 



(5.5) 
(5.6) 

(2) All the minimizing sequences {^n} for E(N) are precompact in H^ s ((M. d ) N ) up to 
translations. Hence there exists {v n } C M. d such that $„(■ — v„) converges, up to a 
subsequence, to a minimizer for E(N) . 



As we will explain in the proof, the role of the additional condition ()5.6[) is to avoid 
vanishing. 



Proof. We split the proof in several steps. We start by proving that the inequalities (|5.5|) 
and (|5.6[) always hold true when the strict inequality < is replaced by a large inequality ^, 
and that if there is equality, then there exists a minimizing sequence which is non-compact, 
for any translations. This shows that (2) implies (1). 



Step 1: Large binding inequalities. The inequalities in (|5.5j) always hold true when the 
strict inequality < is replaced by sj. If there is equality for some 1 ^ A; ^ N — 1, then there 
exists a minimizing sequence {^ n } for E(N) which is not compact, even up to translations. 

Proof. The proof proceeds as usual by constructing a trial sequence ^> n = ^ o ^^(- — R n v) 
(with o = A for fermions and o = V for bosons), where ^ and "J^ are minimizing sequences 
of compact support for E(N — k) and E(k) and R n is large enough. The energy is decoupled 
by (A3). We omit the details. □ 

Step 2: Large inequality (|5.6|) . The inequality (|5.6|) always holds true when the strict 
inequality < is replaced by ^ . If there is equality, then there exists a minimizing sequence 
V&n} for E(N) which is not compact, even up to translations. 

Proof. Removing the center of mass by performing the change of variables x' = 

x\, we see that the original Hamiltonian 



H(N) can be rewritten as 

/ jV-l 



x 2 - Xi, 



N — l 



X N 



H(N) 



bo 



/ 12 



2N 



2N 



E 

7=1 



w 2 



JV-l 

j=i 



JV-l 



E W W 

3=1 



E W(x' k -x' e ) 



H'(N - 1). 
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This shows that the bottom of the spectrum of H(N) is also the bottom of the spectrum 
of H'(N — 1). To account for the original statistics of our particles, the latter Hamilto- 
nian H'(N — 1) is restricted to (N — l)-body functions $ that are symmetric (bosons) or 
antisymmetric (fermions), and additionally satisfy the following relation 



§(-si,a/ 2 



x\) =T<S>(X' 1) X' 2 , 



' N- 



-l) 



with t — \ for bosons and r = — 1 for fermions. Let {$ n } be a Weyl sequence for the bottom 
of the spectrum of the Hamiltonian H'(N — 1) (under the appropriate symmetry constraints) 
and let <p n := n~ d / 2 tp(-/n) for a fixed normalized function ip g H 2 (R d ) n L°°(R d ). We take 
as test function the product state 



whose density is 
P*n {x) — N / dx2 

JR d 

This proves that 



,Xn) = <Pn 



E 



N 



§ n {x 2 • ,X N - Xi) 



dx 



N 




|P*nL 



Id) 



A|M| 



|$„(x 2 - x, ■ 



0. 



,x N 



(5.7) 



hence that ,o^ n —> in L p (M d ) for all 1 < p ^ oo. Under our assumption (Al) on the 
nonlinearity F, this implies that F(p^ n ) 0. On the other hand we have by construction 

lim {V n ,H(N)* n )=wI<r(H(N)) 

71— >oo 

and it follows that E(N) ^ inf a(H(N)). If there is equality, the previous sequence {^ n } 
furnishes a vanishing minimizing sequence. It is not compact, even up to translations. This 
ends the proof of Step 2. □ 



The previous steps show that (2) implies (1). We now turn to the proof of the converse 
implication. We consider a minimizing sequence {^n} an d note that it is necessarily bounded 
in Hl /s ((R d ) N ), by pT3]l . As usual we denote by T n = © • • • © |*„)(*„| the associated 
mixed state in the truncated Fock space. We define like in the proof of Theorem 14.21 the 
number 

m({r„}) := sup jlV(A"r) : 3{v k } C R d , r^r n ,r_^ -+ g A. (5.8) 
We start by proving that vanishing does not hold, that is m({T n }) > 0. 
Step 3: Absence of vanishing. One has m({T n }) > 0. 

Proof. As we have already seen in Lemma l4~71 m({r„}) = is equivalent to having p^ n —> 
strongly in L p ((R™) ) , for all 1 < p < p* . By Assumption (Al), the function F is uniformly 
continuous on L Pl (R d ) n L P2 (R d ) for some 1 < Pl < p 2 < p* . Hence m({r„}) = implies 
that F(p^ n ) and therefore that 

E(N) = lim = lim (V , H (N)V) > inf a(H(N)). 



16 



Mathieu LEWIN 



This contradicts (1), hence shows that it must hold m({r„}) > 0. □ 

Up to a translation (we use that E is translation-invariant) and extraction of a subse- 
quence, we may therefore assume that T n T geometrically, with Tr(A/T) > 0, that is 
r = Goo © • • • © Gnn with ^ Goo < 1- In order to show that {r„} is compact, we have to 
prove that Tr(G NN ) = 1. This only shows that >3>„ strongly in Ll /s ((R d ) N ) but strong 
convergence in H^ s ((R d ) N ) follows by usual arguments. 

Step 4: Decoupling via localization. In this step we split r„ into a part which converges 
to r strongly and a part which escapes to infinity. Contrary to the previous sections, we 
use a radius of localization which depends on n, following Lions [38] [39]. The following is a 
well-known result: 

Lemma 5.1 (Dichotomy). Up to extraction of a subsequence, it holds 



n— voo 



lim / p^ n (x)dx = I pr(x)dx, 



lim / { p^ n {x) + \\J ^/pl^J~x)\ 2 ) dx 

= lim / dxi I dx 2 --- dx N {Wx^^X!, ...,x N )\ 2 = 

n ^°° J R, l ^\x 1 \^6R n JR d JR d 

for a sequence R n — > oo. 

The proof of this lemma uses concentration functions in the spirit of Lions [38] 139] as 
well as the strong local compactness of p*„. See for instance Lemma 3.1 in [TO] for a similar 
result. Let x be a smooth radial localization function with ^ x ^ 1, x( x ) = 1 if |sc| ^ 1 and 
x(x) = if |x| ^ 2, and let 77 := \/l — x 2 ■ Let us consider the smooth localization functions 
Xn '■= x(' I ' Rn) and r/ n — r/(-/R n ), in and outside the ball of radius R n . By Lemma |3. 11 wc 
have (r„) Xii — *- g r geometrically. However by Lemma 15. II it holds 

lim Tr [(L„) X J« = lim / ( X „)V„ = / Pr = Tr [r]«. 

n^oo n^oo J Rd J Rd 

This shows that [(r n ) x jW — » [T]^ strongly in the trace class, hence by Lemma |2~31 that 

(Fn)xn L strongly in S(J^ N ) as n — > 00. 

We can now show that the energy decouples. For the linear part we have by the IMS 
formula (like in the proof of Theorem 13. ip 

CN 

(V n ,H{N)* n ) > Ti>«« (H (r„) x J + Tr^« (H(F„) t) J - — 

+ N(N-1) dan--- / dx N W(x 1 -x 2 )xn(x 1 ) 2 i ln (x2) 2 \^n(xi,...,x N )\ 2 , (5.9) 

JR d JR d 

where H = (B^Li -^( n ) is the second quantization of H{N) in . Performing a 
decomposition similar to (13. lip and using Lemma [5. 11 one sees that the last term of 
goes to zero as n — > oo. For the nonlinear term, we write 

P*„ = |X«| 2 P*„ + |?7«| 2 P*„ = |X«| 2 P*„ + l??n| 2 |X3flJ 2 p*„ + l%fl„| 2 p*„- 
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By Lemma Owe have that |?7 n | 2 |X3i?J 2 P*„ -> in L 1 ^) n L p "(R d ), hence in L Pl (R d ) n 
L P2 (M. d ). Using that F is locally uniformly continuous on L Pl (R d )nL P2 (R d ), we deduce since 
p 9n is bounded in L Pl (R d ) n L P2 (R d ), that 

^(P*J = ^ (l*n| W + l%K„| W) + 0(1). 

By Assumption (A3) we have 

F (IXn| W + l^fljVJ = F (|X»| + F (\ mRn | 2 p*J + 0(1). 

Using again that |??n| 2 |x3-R„ | 2 P*„ — > we finally deduce that 

*Xp*J = F (|x„lW) + F (|77»| 2 P*„) + o(l). 
Hence we arrive at the following estimate 

(tf „, H(N)* n ) > Tr^« (H (r„) x „ ) + F (p(r n)xn ) + Tr^» (H (T n ) Vn ) + F (p (r „ )tta ) + o(l). 

(5.10) 

Let us write the localized states on F^ as 

(r n ) x?1 = Gq' © • • • © G^ , {^n)ri„ = Gq' © • • • © G^ . 

By the concavity of F, we have 

JV 

F(p (rJ J^Ii(Gf)F(pfi r ), (5.H) 
j=o 

with GJ'" := Gj'™/Tr(Gj' n ) (and an obvious convention when G^' n = 0). Using the funda- 
mental relation Tr(G*'™) = Tr(G^" .), we arrive at the lower bound 

N N 

Tr^» (H(r„) r? J + F(p (r , i)rj J > ^Tr(Gf") £(G^V) > ^Tr(Gf l ) F(iV-j). (5.12) 

In the previous bounds, the energy £ is extended to mixed states of Sj N in an obvious 
fashion. Furthermore, for any mixed state G € S($j N ), we have, writing G = ^ Sjl^jX^j I 
with J2j9j = 1 , 

£(G) = 9i{*i,H(N)9j) + W E UP** ) > E 9j > F(N), 

3 \ 3 J 3 

by the concavity of F . Therefore minimizing over mixed states is the same as minimizing 
over pure states, a property that we have used in (|5.12l) . 

Coming back to the term involving Xn m (|5.10p , we claim that it holds 

liminf (ti ><n (H (T„) Xb ) + F (p (r „) x J ) > Tl ><« (HT) + F (p r ) . 

Indeed the interaction term and F(pr n ) converge as n — ¥ oo, by the strong convergence of 
(r n ) x „ towards T in ei(J^ N ). The kinetic energy is lower semi-continuous, by Lemma |2~H 
Summarizing, we have obtained the following lower bound 

N 

E(N) > Tr^N (Hr) + F ( Pr ) +E Tr ( G w) E ( N - ■?')• ( 5 - 13 ) 
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Using the concavity of F as for (T n ) nn , we have 

JV 

Tr^» (WT) + F ( Pr ) > £ Tr(G w ) E(j), 
j=o 

hence it follows that 

N 

e(n) > £ Tr ( G «) ( E ti) + E ( N - i))- 

When the binding condition ()5.5[) holds true, this is only possible when Gu = • • • = 
Gn-i n-i = 0. 

Step 5: Conclusion. It rests to prove that Goo = 0. Let * be the weak limit in $) N of the 
original minimizing sequence {*f> n } and notice that Gat at = |$)($|. Since Gat at ^ 0, it holds 
W 7^ 0. Inserting all this in (|5.13l) (recall pa 00 — 0), we obtain the estimate 

(1 - Tr(G oo ))£(A0 = \\yfE{N) > (% H(N)^} + F(p*). (5.14) 

If 1^1 < 1, then we use (A5) and get 

F(p») > ||*|| 2 F(^ /W |), 

that is 

(%H{N)*) + F(ft V ) > ||*|| 2 £ (j^j > ll^ll 2 ^^). 

This contradicts (|5.14p . hence implies that it must hold ||$[ = 1 and Goo = 0- This ends the 
proof of Theorem 15.11 □ 

Theorem 15.11 can be generalized to finite-rank fermionic systems (Hartree-Fock case or 
multiconfiguration theory when W 0), following the arguments of Section[4] For instance, 
in the Hartree-Fock case one can easily prove the following 

Theorem 5.2 (Nonlinear HVZ for many-body systems in the Hartree-Fock ap- 
proximation). Let En{N) be the (fermionic) ground state energy in the Hartree-Fock ap- 
proximation, defined by 

E N (N):= inf £(*). (5.15) 

rank(*)=A r 
1*1=1 

Under the previous assumptions, the following assertions are equivalent: 

(1) One has 

E N (N)<E N - k (N-k) + E k (k) for all k = 1, ...,N - 1; (5.16) 

(2) All the Hartree-Fock minimizing sequences {^n} for Em(N) are precompact in 
H^((M. d ) N ) up to translations. Hence there exists {v n } C M. d such that *„(• — v n ) con- 
verges, up to a subsequence, to a minimizer for En(N). 

Note the absence of a condition of the form (|5.6p : as we have seen in the proof of Theorem 
14.21 in the case of vanishing of a Hartree-Fock state, the interaction energy always tends to 
zero. The condition (|5.16p is sufficient to avoid this. 
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5.2. Application: the multi-polaron 

In this section we study a system of N electrons in a polar (ionic) crystal, called N -polaron. 
Thanks to the underlying deformations of the crystal, the N electrons can overcome their 
Coulomb repulsion and form a bound system. Recently there has been a renewed interest in 
the multi-polaron problem, triggered by the possibility of bipolaronic superconductivity in 
high-temperature superconductors |15j . 

Under the assumption that the polaron extends over a region much bigger than the typical 
spacing between the ions of the crystal, one can use a continuous model based on phonons. A 
model of this form was proposed by H. Frohlich in [21] , It assumes a linear coupling between 
the electrons and the longitudinal optical phonons, together with a constant dispersion 
relation for the phonons. The corresponding Hamiltonian takes the form 

E ( Z ^T ~ V^i)) + E I ~ I + / dkaHk)a(k), (5.17) 



where 



2tt J r3 \k\ 

The Hamiltonian acts on the Hilbert space L 2 a ((R 3 ) N ) (g> F s , with a!{k) and a(k) being the 
creation and annihilation operators (in the Fourier representation) for the phonons on the 
bosonic Fock space J- s . Because of its relation to the dielectric constants of the polar crystal 
[20] 164) . the parameter a must satisfy the constraint a < U in the physical regime. For 
simplicity we have discarded the spin of the electrons. 

In the regime of strong coupling, the model reduces to the so-called Pekar-Tomasevich 
(PK) theory [47l 08J [44] in which the interaction with the crystal is modelled by a classical 
Coulomb self-interaction. The energy is now given by 

(5.18) 

for $ e ^^((R 3 )^). The corresponding ground state energy is as usual defined as 

E aiU {N)= mf £ a ,u(?). (5.19) 
11*11=1 

We have emphasized the dependence in the parameters a and U. A simple scaling argument 
shows that E a jj = V E a /j/ t i, hence we may work in a system of units such that U = 1. In 
this case, for simplicity we use the notation £ a := E a .\ and E a (N) :— E aj \(N). 

Another way to derive the Pekar-Tomasevich energy is to restrict to (uncorrelated) prod- 
ucts states of the form $ ® $ g L^((U. 3 ) N ) ® J" s and to minimize with respect to the state 
$ of the phonons [22]. 

Both the original model of Frohlich and the Pekar-Tomasevich theory have stimulated 
many works. On the mathematical side, the validity of PK theory in the large coupling regime 
was shown for TV = 1 by Donsker and Varadhan in [14] . and with a different approach by 
Lieb and Thomas in [37]. The case N = 2 was treated by Miyao and Spohn in [44]. The 
stability or instability of large polaron systems was studied by Griesemer and M0ller [22] , 
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then by Frank, Lieb, Seiringer and Thomas [TBI [17]. In this latter work, the absence of 
binding of iV-polaron for small a is also proven. 

Using geometric techniques, we are able to study the existence of multi-polaron systems: 

Theorem 5.3 (Binding of Pekar-Tomasevich multi-polarons). Assume U = 1. For 
every N ^ 2, there exists a constant t c (N) < 1 such that the following hold for alia > t c (N): 

(1) E a {N) < E a (N - k) + E a (k) for all k = 1,...,N-1; 

(2) All the minimizing sequences {^n} for E a (N) are precompact in H^ s ((BL d ) N ) up to 
translations. Hence there exists {v n } C M d such that ^> n (- — v„) converges, up to a 
subsequence, to a minimizer VP for E a (N); 

(3) Any such minimizer satisfies the following nonlinear eigenvalue equation: 



where fi is the first eigenvalue of the many-body Schrddinger operator in the parenthesis. 

Our result covers the physical range a £ (r c (JV), 1) but we do not provide any bound on 
the critical t c (N). It was proved in [T5] that binding does not occur when a is small enough, 
hence one must have t c (N) > 0. We expect that t c (N) — > 1 when N — > oo but we do not 
have a proof of this. 

For N = 1, the Pekar-Tomasevich energy is defined as 



and it is sometimes also called the Choquard functional. The existence and uniqueness of a 
ground state up to translations for all a > was proved by Lieb in [31] . Nothing seems to 
be known on the uniqueness of ground states up to translations for N ^ 2. 
For the bipolaron (N = 2), the binding energy 



is a convex and non-decreasing function of a. We deduce from Theorem 15 .31 that there exists 
r c (2) < 1 such that binding does not hold for all ^ a ^ t c (2), whereas binding holds true 
and minimizers exist for all a > t c (2). A result of the same form was already announced in 
[44) . Numerical computations 65, 59 suggest that, for the bipolaron, r c (2) ~ 0.87. 

Since the Pekar-Tomasevich model is exact in the limit of strong coupling, a/U < 1 and 
a>l, our result implies the existence of binding for Frohlich's iV-polaron described by the 
Hamiltonian (|5.17|) . when r c (N) < a/U < 1 and a is large enough. For small a, numerical 
computations indeed suggest that Frohlich's polaron does not bind for any U > a. In [65] 
(Fig. 4) the critical value above which Frohlich's bipolaron formation is possible was found 
to be a ~ 13.15. 

Remark 5.1 (Extensions). For anisotropic materials, one can take F of the form 



where M is a 3 x 3 symmetric matrix satisfying M ^ 1 . Existence of ground states follows 
from our method when M is sufficiently close to the identity matrix. 




(5.20) 




(5.21) 



2E a {l) - E a {2) = 2E 1 {1) a 2 - E a {2) 




Geometric methods for nonlinear many-body quantum systems 



51 



Our results hold the same in 2D, assuming the particles interact with the 3D Coulomb 
potential, a model which is often considered in the physical literature (see, e.g. \65l \ 64V ). o 

Thanks to Theorem 15.11 the proof of Theorem 15.31 is essentially reduced to showing the 
binding condition. This is done by building suitable trial states. The easy case is a > 1, when 
two multi-polaron always have a Coulomb attraction at large distances. The case a = 1 is 
more subtle, and we prove that there is always a Van Der Waals attraction at large distances, 
following Lieb and Thirring [36 . The existence of t c (N) is then obtained by continuity of 
a i y E a (N), using that there are only finitely many binding conditions to verify. 

Proof. The energy £ a is of the general form which we have considered in Section I5TT1 The 
nonlinear functional 

F(P) = ~[ I P ~^d X dy 
2 Jr3 Jr3 \x - y\ 

is clearly strictly concave, and it satisfies our assumptions (A1)-(A5) with pi — P2 — 6/5, 
by the Hardy-Littlewood-Sobolev inequality |34) . Furthermore, the condition (I5.6[) reduces 
to E a (N) < since the interaction potential W(x) — l/\x\ is non-negative. This condition 
is implied by the binding condition, hence it is only necessary to verify that E a (N) < 
E a {N — k) + E a (k) for k = 1, N — 1. Since the function a i— > E a (N) is clearly continuous, 
it is sufficient to show that 

E a (N) < E a (N -k) + E a (k) for all integers 1 < k < N - 1 and all a > 1. (5.22) 

As usual, we prove these binding inequalities by induction, assuming that E a (k) has a 
minimizer for all k = 1, N — 1. For N = 1, we already know that ground states of E a (l) 
exist for all a > 0. The following will be very useful. 

Lemma 5.2 (Properties of multi-polaron ground states). Assume that ^ is a ground 
state for E a (N) with a > 0. Then VP solves the self- consistent equation (|5.20j) where fj, is 
the first eigenvalue of the many-body operator 

If a > 1 — 1/N, then fi < inf a css (H^(N)) and both ^ and V^P decay exponentially at 
infinity. 

Proof. We have already explained in the proof of Theorem 15 . 1 1 that . by the concavity of F, 
E a (N) is also the lowest energy over all mixed states of ^((R 3 )^). In particular it holds 

£j(l-t)\y)(y\+t\V)(y'\ \ > E a (N) 

for all e Hl((R 3 ) N ) and all ^ t ^ 1. The first order in t provides the bound 
(*', fff (N)W) ^ H$(N)V}, showing that n is the first eigenvalue of H$(N). The Hamil- 
tonian H$(N) is a usual Coulomb Hamiltonian of N electrons with an external Coulomb 
field of total charge Z = a L 3 = aN. It was shown by Zhislin and Sigalov [66l [67] that 
fi is an isolated eigenvalue as soon as N < Z + 1 = aN + 1. The exponential decay follows 
from well-known results reviewed for instance in Section XIII. 11 of 1521 . □ 
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Let us now assume that E a (N — k) and E a {k) have respective ground states V^i and 
^2, and that a > 1. We want to prove that E a (N) < E a (N - k) + E a (k). Using their 
exponential decay, we can replace and ^ by functions with support in a ball of radius 
R, making an error in the energy of the form e~ aR . For the sake of simplicity we do not 
change our notation and assume that 

£a(*i) < E a (N - k) + Ce- aR , £ a (tf 2 ) < E a {k) + Ce- aR . 

When a > 1, we can take advantage of a Coulomb attraction at infinity and choose as 
trial function 

for rotations U, V £ 5*0(3) and with rotated ground states ^ = ^(C/" -1 -). Averaging over 
the rotations U, V £ 50(3) and using Newton's theorem yields a bound 



/ dU [ dV£ a ^ U R V )^E a (N-k) + E a (k)- {N fc f p (1 a) +Ce-°i< 

JSO(3) JSO(3) 6K 



This shows the binding inequality when a > 1. 

When a = 1 there is a priori no simple binding in l/R. Fortunately, there always exists a 
Van Der Waals force between two multi-polarons. Following a method of Lieb and Thirring 
|36) . we take as trial state 

:= *i A - 3ifi7) + A < | m ■ > ; V; | } A { | n ■ > ; V, | ( 3/?D 



Writing with an obvious convention — + ^Y > we have 




dx 2 ■ ■ ■ / dx N ^Y^ U r V = 

which is seen by using that "J/f and ^ v (■ — 3i?u) have disjoint supports, as well as the fact 
that is orthogonal to ■ J^^i' ^i) ^i 7 anc ^ a similar property for ty v ■ As was already 
mentioned in [36], this yields ||^jj V '|| 2 = 1 + 0(A 2 ), but this also gives 

p K ,v /\\*r v \\ ^ +/ 3 * 2 v (-- 3 - r «) + °( a2 )- ( 5 - 23 ) 

Therefore we can mimic the argument of [36] and obtain an upper bound of the form 
/ dU I dV £ a fa% v '/ )^Ex{N-k) + E x {k) + a-^ + b\ 2 + C e - aR . 

JSO(3) JSO(3) V / R 

The linear term in A comes from the cross-term between the two functions appearing in the 
definition of vf^' , in the electron-electron interaction term. This term is exactly the same as 
the one calculated in |36j . The nonlinear term involving the density only provides a 0(A 2 ) 
by (|5.23|) . Taking A = — a/26i? 3 yields the desired attractive Van Der Waals interaction 
potential — O/i? 6 , hence the binding of two polaron systems when a = 1. This ends the 
proof of Theorem 15.31 □ 
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